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List  of  fiyrobolB  '  ' 


r,  8,  z  Cylindrical  coordinates,  8ee  Fig*  (l-l) 

ufz,  0,  t),  v(z,  8,  t)  Longitudinal,  tangential  and  radial  cctt- 

v(t,  8,  t)  ponents  of  the  shell  displacement*  (v  is 

measured  positive  inward),  Bee  Fig. 


u(t,  0),  v(z,  0),  v(z,  0) 

a 

Vj*  J  *  0-8 


Space  dependent  parts  of  u,  v,  and  V  Is  S 
principal  mode. 

Radius  of  shell 

Fourier  Series  coefficient  -  Bee  Eq.  (VII-18)* 

(VII-21) 


Generalized  coordinate  -  sloshing  modes  of 
fluid  filling  in  shell. 


Cross  sectional  area  of  wind  girder. 
Normalization  Coefficient 

Expansion  coefficient  for  velocity  potential 
Coefficient  -  8ee  Eq.  (lV-52)-(lY~53). 
Constants  -  See  Eq.  (l-l)-(l-3)» 

Modulus  of  Elasticity 
Integral,  See  Eq.  (VII-15). 

Coefficient  -  See  Eq.  (IV-95)* 

Shell  Thickness 
Integral,  8ee  Eq.  (VII-32) 

Moment  of  Inertia  of  wind  girder. 

BeBsel  Function  of  first  kind  of  order  n. 

Coefficient. 


Mass  ratio,  K 


mvn 


Vs 


(l)  Additional  symbols  in  the  text  are  defined  as  they  occur. 
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r, 
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M, 


M_ 


Mnl 
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1 


m 


vn 


n 


P 

P(e,  *,  t) 
P(r,  e,  t) 


^(t) 

V  Sol*  Si2 

8k 

% 

t 

T 

T 


Height  of  shall  . 

Height  of  wind  girder  shove  bottom  of  tank. 

Frequency  number.  See  Eq.  (1-22). 

Generalized  mass  of  shell  in  mode  n  »*  0, 

See  Eq.  (IV-79)* 

Total  mass  of  empty  tank  ■  2«alm^  . 

Frequency  numbers  -  Sea  Section  (ill). 

Mass  per  unit  of  area  of  the  empty  shell. 

th 

Virtual  maps  per  unit  ares  of  the  fluid  in  the  n 
mode,  See  Eq.  (IV-55) * 

Number  of  circumferential  waves  of  a  mode  of 
vibration. 

Pressure  in  fluid. 

Radial  blast  pressure  on  shell. 

Vertical  blast  pressure  on  fluid  surface 
or  on  roof  of  shell 

Radial  blast  pressure,  in  subarea  "ih".  See  Section 

(Vll-a) 

Vertical  blast  pressure  on  subarea  8^. 

Generali zed  coordinate  of  shell  displacements. 

Generalised  forces  in  the  nth  mode. 

Sector  of  shell,  See  Section  (VIX-a). 

Subarea  over  which  vertical  blast  forces  on 
the  fluid  surface  of  the  shell  are  measured. 

Time 

Kinetic  energy 
8ee  Eq.  (V-39) 

Kinetic  energy  of  fluid 
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Coefficients  defining  the  shapes  of  modes 
of  free  vibration  of  the  shell. 

Potential  Energy  of  Shell. 

Fluid  velocities  in  the  r,  9  and  z  directions. 

Root  of  transcendental  Eq.  (IV-23). 

Coefficient  -  See  Eq.  (IV-62)  and  (IV-93). 

Coefficient  -  See  Eq.  (V-112) 

Conversion  Coefficient  -  See  Eq.  (V-42)-(V-45) 

Fraction  denoting  the  height  of  the  fluid 
filling  in  the  shell.  Height  of  fluid  =  yL. 

Coefficient  -  See  Eq.  (IV-58). 

Shell  Strains. 

Displacement  of  Fluid  Surface. 

Poisson’s  Ratio. 

Ratio  of  height  of  tank  to  radius  of  tank. 
Mass  density  of  fluid. 

Mass  density  of  empty  shell. 

Shell  Stresses. 

Velocity  potential  function  of  fluid 
See  Eq.  (VIX-33). 

Frequency  of  vibration  of  empty  or  partially 
full  shell. 


Note:  r,  9,  or  z  used  as  a  subscript  for  the  shell  displacements  u,  v, 
v,  denotes  differentiation  with  respect  to  the  particular  variable  used. 


Dots  indicate  differentiation  with  respect  to  time 


Consider  the  empty  fuel  tank  of  1?ig.  (i-i)  consisting  of  a  thin 
cylindrical  shell  of  height  L  and  radius  "a')  free  at  the  top  and  simply 
supported  at  the  base.  The  positive  directions  of  the  displacement 
components  u,  v  and  w  are  3hown  in  the  figure. 

An  analysis  of  the  frequencies  and  modes  of  such  a  structure  has 
been  made  in  Reference  (i)^  on  the  assumption  that  the  shell  acts  as  a 
membrane  with  no  bonding  stiffness.  In  addition,  a  procedure  for  obtain¬ 
ing  the  frequencies  of  the  shell  including  bending  effects  is  presented 
in  this  reference.  The  condition  of  simple  support  at  the  base  of  the 
shell  requires  in  general  that  the  displacements  u(z,0,t),  v(z-0,t)  and  w(z ,0 

t)  and  the  moment  M  be  zero  at  z  =  0.  If  the  shell  is  considered  to  be 

z 

a  membrane  with  no  bending  stiffness,  only  the  conditions  that  u(z,  0,  t) 

2 

and  v(z,  0,  t)  are  zero  at  z  =  0  can  be  enforced. 

The  transcendental  equations  for  the  frequencies  of  vibration  and 
the  expressions  for  the  mode  shapes  corresponding  to  these  frequencies, 
according  to  Reference  (l),  are  extremely  complicated  and  ore  very  diffi¬ 
cult  to  apply  to  forced  vibration  problems  involving  the  response  of  the 

(1)  "Dynamic  Response  of  Cylindrical  Tanks"  by  P.  L.  DiMaggio,  Armed  Forces 
Special  Weapons  Project,  Contract  DA-2y-004~XZ-54,  AFSWP  No.  1075,  May,  1958. 

(2)  Any  combination  of  two  displacements  set  equal  to  zero  could  suffice 
to  solve  the  boundary  value  problem  for  the  membrane  shell.  However,  as 
shown  in  Reference  (l),  the  particular  combination  used  gives  rise  to  very 
small  displacements  w(z,  0,  t)  at  the  base  of  the  shell.  The  correction 

to  the  displacement  curve  for  w(z,  0,  t)  when  bonding  stiffness  is  included 

in  the  analysis  and  the  condition  w(z,  0,  t)  =  0  at  z  =  0  is  enforced  is 
shown  to  be  very  small  under  these  conditions. 
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shell  to  transient  blast  loadings.  In  the  present  section  of  this  report, 
using  these  frequencies  and  shapes  as  a  guide,  approximate  values  for  the 
frequencies  and  modes  of  free  vibrations  of  thin  cylindrical  shells  of 
constant  thickness  are  obtained.  These  approximate  modes  are  of  a 
relatively  simple  form  and  can  easily  be  used  in  Section  VII  of  this  Report 
in  which  an  analysis  of  the  forced  vibration  of  empty  and  partially  full 
tanks  is  presented.  Since  the  application  of  these  results  to  both  steal 
shells  and  steel  shells  with  concrete  shielding  is  contemplated,  both 
membrane  and  bending  effects  will  be  considered  in  the  following  analysis. 

The  modes  of  free  vibrations  of  the  shell  may  be  characterized  by 
an  integral  number  n^  the  number  of  circumferential  waves  in  the  mode. 

For  each  configuration  there  exists  an  infinite  number  of  frequen¬ 
cies  and  corresponding  mode  shapes.  In  general,  for  n  /  0,  we  will  only 
require  the  lowest  one  of  these  frequencies  and  mode  shapes.  For  the 
case  n  •  0,  special  conditions  prevail,  and  this  case  will  be  discussed 
in  Part  (  e  )  of  this  section. 

The  displacements  u,  v,  w  of  the  shell  in  free  vibrations  can  be 
expanded  into  a  power  series  in  the  coordinate  zt^2a^ 


u(z,  9,  t) 


_  z  _  z  „  z 

C  —  +  C  — ~  +  C  — ~  +  .  •  < 

u,  a  uc  . 2  uQ  a3 

1  2  a  3  a 


coanS  e 


lust 


(I-D 


V(t,  9,  t)  a 


Cv  a  +  Cv  ^2  +  Cv  \  +  ••• 

V1  a  v2  a2  v3  a3 


uinnfi  ei<Dt  (l»2) 


v(z,  9,  t)  - 


Cw  +  c_.  i  +  c 


a 


v  2 
2  a 


cosnd  e 


icot 


(1-3) 


(2a)  From  the  reference  of  Footnote  (l)of  this  section,  it  is  evident  that 
the  z  dependent  terms  of  the  mode  shapes  are  well  behaved  functions  which 
can  be  expanded  into  power  series  about  the  point  z  *  0. 


.  ,*  -\  «\  •*.  -\  -v  •> 


where  the  constant  terms  C  and  C  are  taken  equal  to  zero  due  to  the 

uo  vo 

boundary  conditions  at  the  base  of  the  shell.  In  the  membrane  analysis 

as  presented  in  Reference  (l),  the  radial  displacement  cannot  be  zero 

at  z  »  0  and  the  constant  term  C  must  be  included  in  the  series  for 

o 

w(z,  6,  t).  The  choice  of  the  number  of  terms  in  each  series  which  must 
be  retained  in  an  spproximate  analysis  is  governed  by  the  simplest  combina¬ 
tion  of  terms  required  to  obtain  good  approximations  to  the  frequencies, 
membrane  strains  snd  membrane  stresses  of  Reference  (l).  It  should  be 
noted  that  the  exact  membrane  displacements  do  not  differ  very  much  from 
straight  lines  and  hence,  relatively  few  terms  need  be  retained  in 
8q«  (1-1)  •  Xq.  (1-3).  A  detailed  discussion  of  the  above  considerations 
appears  in  Section  II  of  this  Report. 

Two  seta  of  approximate  displacements  are  considered  in  this  section; 

1)  a  set  in  which  five  constants  are  retained  in  Iq.  (I-l)  to  (1-3)  end 

2)  a  set  in  which  three  constants  are  retained.  While  the  latter  will 
give  satisfactory  estimates  of  frequencies  for  most  applications,  the 
former  are  required  to  give  accurate  values  of  strains  and  stresses  for 
use  in  forced  vibration  problems. 


where  n  *  1,  2,  ....  is  on  integral  number  indicating  the  number  of 

circumferential  waves  in  the  mode,  <o  is  the  frequency  of  free  vibrations 

and  U,  V,  W,  X  and  Y  are  the  five  constants. 

(2c) 

The  Royleigh-Ritz  method  will  be  used  to  determine  the  frequencies 
and  mode  shapes  of  the  free  vibrations  of  the  shell.  To  use  this  method, 
the  kinetic  and  potential  energies  of  the  shell  must  be  determined  in 
terms  of  the  parameters  U,  V,  W,  X  and  Y.  The  kinetic  energy  of  the  shell 
is  given  by  the  relation 

2n  L 

T  *  sT  J J  +  +  ^  ad0ds  (1-7) 

o  o 

where  is  the  maBB  per  unit  area  of  the  shell 

mi  *  ^i*1  (1-8) 

and  dots  indicate  differentiation  with  respect  to  time. 

(2b)  The  expressions  in  the  brackets  in  Eq.  (1-4)  and  (1-6)  for  u  and 

p 

w  are  of  the  form  (az  +  bz  )  and  (c  +  dz)  respectively.  The  particular 


(2b.  Cont.).  .  .combinations  of  coefficients  used  in  these  equations 
were  chosen  so  that  the  kinetic  energy  of  the  empty  shell,  Eq.  (I —9) > 
would  contain  only  the  squares  of  the  coefficients  U,  V,  W,  X  and  Y. 

This  leads  to  a  considerable  simplification  in  the  frequency  determin¬ 
ant,  Eq.  (1-23) ,  by  making  the  YW  and  WY  coefficients  equal  to  zero  and 
insuring  that  all  terms  containing  the  frequency  number  M  lie  only  on 
the  main  diagonal. 

(2c)  The  Rayleigh-Ritz  Method  is  based  on  the  condition  that  the  varia¬ 
tion  of  the  difference 

L  -  T-V 


vanishes  for  arbitrary  variations  of  the  deflection  curve  consistent 
with  the  boundary  conditions,  i.e.  the  geometrical  restraints,  of  the 
problem.  By  chosing  the  displacements  to  be  functions  containing  the 
arbitrary  coefficients  U,  V,  W,  X  and  Y,  the  quantities  T  and  V  can  be 
expressed  in  terms  of  these  coefficients,  See  Eq.  (l-9)-(I-l4)  and  the 
function  T-V,  may  be  formed.  This  quantity  is  now  a  function  of  the 
variables  U,  V,  W,  X,  and  Y.  The  condition  that  the  variation  of  T-V 
vanishes  for  arbitrary  variation  in  U,  V,  W,  X  and  Y  requires  that 


where  takes  on  successively  the  values  U,  V,  W,  X  and  Y.  This  leads 
to  a  system  of  five  simultaneous  homogeneous  equations  on  U,  V,  W,  X  and 
Y;  for  non  zero  solutions  of  this  system  of  equations,  the  determinant 
of  the  coefficients  must  be  set  equal  to  zero  thus  leading  to  the  frequency 
determinants  which  appear  later  in  the  paper.  See,  for  example,  "Vibration 
Problems  in  Engineering"  by  S.  Timoshenko,  D.  VanNostrand  and  Co.,  Third  Edition, 

Jan.  1955,  Pg.  381  ff  •  In  this  reference  the  method  is  called  just  the  Ritz 
method.  Also  "Mathematical  Methods  in  Engineering"  by  T.Von  Karman  and 
M.  Biot,  McGraw  Hill,  19*t0,  Pg.  352  ff. 
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(3)  "Tables  for  the  Frequencies  and  Modes  of  Free  Vibration  of  Infinitely 
Long  Thin  Cylindrical  Shells"  by  M.  L.  Baron  and  H.  H.  Bleich,  Journal  of 
Applied  Mechanics,  Vol.  21,  No.  2,  June  1954.  Transactions  of  the  American 
Society  of  Mechanical  Engineers. 


l7~)  7  J oi  L  “ +  M  +  v)  +  <¥>  (“ 8 '  “8) 


+  21 |=JiI  »2(vt  ♦  „iS)2  +  2a2V»It(wes  +  re 


+  2a- 


w  ( 

t  ti 


The  subscript e  Indicate  partial  differentiation  with  respect  to  %  or  B, 
Substituting  the  expressione  for  u,  v  and  w  into  Iq«  (i-ll)and  Iq.  (1-12) 
the  potential  energy  of  the  shell  is  given  bys 

o2  [i  +  t  v2  [lirJil  +  a^]+  v2^ 

V  „  I&5 _  | 

1  £(1-/)  ^X2  |  t2  ♦  i±^SX|  ♦  Y2  ♦  UV  j^i)  nT] 

+UX  l  -  UY2V  -  VW  -  VYnl 


(4)  It  is  proper  to  call  attention  to  the  fact  that  the  expression  for 
the  strain  energy  Vg  used  may  not  be  rigorously  correct  in  the  light  of 
recent  advanced  shell  theories'*.  However,  the  authors  are  satisfied  that 
the  effeat  of  any  corrected  strain-energy  egression  in  the  range  of  appll 

cation  of  these  results  is  insignificant. 

(5)  "The  New  Approach  to  Shell  Theory",  By  X.  H..Kennard,  Journal  of 
Applied  Mechanics,  Vol.  75*  1953*  Transactions  of  the  American  8ociety 
of  Mechanical  Engineers  ,  pp  33-40. 


where  ft  «  L/a. 

The  Rayleigh-Ritz  method  for  the  determination  of  the  frequencies  o> 
and  the  ratio  of  the  shapes  requires  that 

33  (T  -  \  -  v2)  «  o  (1-15) 

where  successively  takes  the  values  U,  V,  W,  X  and  Y.  Eq.  (1-15)  leads 
to  the  following  set  of  five  homogeneous  linear  equations  in  the  five 
unknowns  U,  V,  W ,  X  and  Y: 
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1  -  k(l-n  ff  Y  -  0  (1-20) 
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Nonvaniahing  solutions  of  Bq.  (1-17) -( 1-20)  and  free  vibrations  exist  only 
if  the  determinant  of  the  system  vanishes.  This  leads  to  the  determlnant&l 
frequency  equation,  shown  on  Page  13. 

For  given  values  of  |  and  n,  JSq.  (1-23)  yields  five  positive  roots 
My  defining  five  mutually  orthogonal  modes,  of  frequencies 
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J  -  1,  2,  3,  4,  5  (1-24) 
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Membrane  +  Bending  Effects 
5  Constant  Approximation 
Empty  Shell 


The  shape  of  the  mode  pertaining  to  a  particular  frequency  a>j  can  be 

U  V  X  Y 

found  by  computing  the  ratios  ^  -  and  -  from  any  four  of  the  Eq. 

(l-l6)-(l-20) .  In  general,  the  lowest  frequency  only,  will  be  required 

and  for  simplicity  in  the  following  expressions  the  subscript  j  for  each 

of  the  constants  U.,  V.,  W  ,  X  and  Y  will  be  dropped. 

J  J  <3  J  J 

Once  the  frequency  a>  and  the  ratios  of  the  constants  have  been  ev¬ 
aluated  for  a  particular  value  of  n,  the  displacements  of  the  shell  may  be 
written  in  the  form: 


u(,'8)  ■  Cn  H  +  S  -  ft >  COBn$ 

a  4a 


v(z,  e) 


C  ^  —  sinnQ 
n  rf  a 


w(z,  0)  “  cn  w  +  (f  -  §“)  cosnS 


It  is  convenient  to  define  "normalized  modes  of  vibration"  by  chosing  the 


constant  C  so  that 
n 

L  2n 


/o/o(uS  +  v8  +  w^)adOdz  =  2jtaLm^  *  M., 


In  thiB  case,  the  modes  are  normalized  to  the  total  mass  of  the  empty 
shell,  M^.  Substituting  Eq.  (l-25)-(l-27)  into  Eq.  (1-20),  the  normaliza¬ 
tion  coefficient  C  becomes 
n 
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f  A2  4 .  /Vx2  ,  i  ,  3/Y\2  3.2,XvJ 

~¥  +  ¥  +i  +  p(w)  +  Bo6  ¥ 


UHing  Eq.  (1-29),  Eq.  (1-25) -(1-27)  give  the  normalized  modes  of  free 
vibration  of  the  empty  tank. 

It  ie  of  interest  to  evaluate  the  normalized  Btraina  and  the  correo- 
ponding  normalized  stresses  in  each  mode.  The  strains  for  the  plane 
stress  solution  of  the  boundary  value  problem  for  the  free  vibrations 
of  thin  cylindrical  shells  are  given  by 
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The  normalized  stresses  may  then  be  obtained  from  the  relations 
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b)  Approximation  Using  Five  Constants  -  Membrane  Effects  Only 
n  £  0 

For  the  case  of  thin  shells  in  which  the  ratio  of  the  thickness  to 
the  radius,  -,  is  very  small,  the  strain  energy  of  bending  V0  becomes 

Cl  Cm 

extremely  small  with  respect  to  the  membrane  strain  energy  V^,  for  modes 

(6) 

with  a  low  circumferential  wave  number  "n'1'  *  In  such  cases,  excellent 

approximations  of  frequencies,  strains  and  stresses  can  be  obtained  if 
the  shell  is  considered  to  act  as  a  membrane  with  no  bending  stiffness 
and  Vg  is  set  equal  to  zero.  The  range  of  h/a  and  n  for  which  the  membrane 
assumption  is  permissible  is  discussed  in  detail  in  Section  III  of  this 
report  for  several  cases  of  interest. 

The  frequency  equation  and  the  mode  shapes  for  the  membrane  shell  are 
obtained  by  setting  the  value  of  the  constant  k  equal  to  zero  in  Eq.  (l-l6)- 
(1-20).  The  system  of  homogeneous  equations  then  becomes: 


(1-39) 


(6)  It  should  be  noted  that  for  higher  modes  with  very  large  values  of  n, 
the  effect  of  the  bending  strain  energy  Vg  becomes  prevalent  and  the  frequency 
in  the  higher  modes  is  controlled  by  the  inextensional  effects. 


»  «  *■  *<  »  w  *  *  V  -  »  -  *  -  * 

/>  >»-•  V.VA&.W  .Vf V ’.V 


The  determinantal  frequency  equation  is  obtained  by  setting  the  determin¬ 
ant  of  the  above  system  equal  to  zero,  or  by  setting  k  ■  0  In  Xq.  (1-23) • 
This  determinant  Is  shovn  on  Page  18. 

For  given  values  of  {  and  n,  Xq.  (X-4l)  yields  the  five  positive 
roots  My  The  shape  of  the  mode  pertaining  to  a  particular  Mj  can  be  found 
by  computing  the  ratios  g  ,  g,  g  and  g  from  any  four  of  the  Xquatlons  (1-36) 
to  (1-ltO).  The  frequency  can  be  computed  from  Iq»  (1-24).  In  general, 
only  the  lovest  frequency  <n^  will  be  required. 

Once  the  frequency  m  and  the  ratios  of  the  constants  have  been  ev¬ 
aluated  for  a  particular  value  of  n,  Xq.  (I-25)-Kq»  (1-35)  may  be  used 
to  compute  the  normalised  mode  shapes,  normalised  strains  and  normalised 


stresses  In  the  shell 


c)  Approximation  Using  Three  Constants  -  Membrane  and  Bending  Effects. 

n  4  0 


Let  the  displacements  of  the  shell  (Fig.  I-l)  he  given  by 
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v(tt  0,  t)  «  v(z,  ©Je*^  a  W“  cos  n6  e 

0L 


icut 


(1-44) 


Substituting  Eq.  (l-42)-(l-44)  into  Eq.  (1-7) >  the  kinetic  energy  of 
the  shell  is  given  by 
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The  potential  energy  of  the  shell  is  again  given  by  Eq.  (l-10)-(l-3L2), 
which  vqpon  substitution  of  Eq.  (1-42) -(1-44)  become: 
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Using  Sq.  (I“15),  the  folloving  system  of  three  homogeneous  equation# 
in  the  three  unknown*  U,  V  and  W  is  obtained: 

[zk  -  4  -  n2a-V)(Xtloj  0  +  V  w  •  0  <I*M) 
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where  k  and  M  are  defined  by  Iq.  (1-21)  and  (l-22)  respectively.  Honvanish- 
ing  solutions  of  the  system  of  Kq.  (1-48) -(1-50)  and  free  vibrations  exist 
only  if  the  determinant  of  the  system  vanishes.  This  leads  to  the  deter¬ 
minants!  frequency  equation  shown  on  Pg.  21. 

for  given  values  of  |  and  n,  Kq.  (IV- 51)  yields  three  positive  roots 
defining  three  mutually  orthogonal  nodes  of  frequencies 

o  M.Kh 

<£  .  -i- -  J  -  1,  2,  3.  (1-52) 

4  a1a2(l-V2) 

The  shape  of  the  mode  pertaining  to  n  particular  frequency  ct^  can  be 
found  by  computing  the  ratios  ^  and  ^  from  any  two  of  Iq.  (l-48)-(l-50). 

In  general,  the  lowest  frequency  cc^  only,  will  be  required  and  for 
simplicity  in  the  folloving  expressions,  the  subscript  J  will  be  dropped. 
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The  displacements  of  the  shell  may  be  written  in  the  form 


u(z,  6)  m  Cn  §  |  cosn0 


(1-53) 


v(s,  e)  m  Cn  ~  ~  sinne 


(1-54) 


w(t,  6)  m  C  —  coend. 
n  a 


(1-55) 


It  is  convenient  to  define  "normalized  modes  of  vibration"  by  chosing 
the  constant  Cn  so  that  Eq.  (1-28)  is  satisfied.  In  this  case,  the  modes 
are  normalized  to  the  total  mass  of  the  empty  shell,  ■  2nalm^ .  Sub¬ 
stituting  Eq.  (l-53)-(l-55)  into  Eq.  (1-28),  the  normalization  coefficient 

C  becomes 
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.1 


,  /U\2  /Vv2 

1+(w)  +<W)-U 


k 


(1-56) 


Using  this  value  of  Cq,  Eq.  (l-53)-(l-55)  give  the  normalized  modes 
of  free  vibration  of  the  empty  tank. 

The  results  obtained  from  this  approximation  may  be  used  where  an  esti¬ 
mate  of  the  frequency  of  vibration  in  any  particular  mode  is  required.  The 
work  required  in  solving  the  third  order  determinantal  frequency  equation 
is  materially  less  than  that  required  for  the  fifth  order  determinantal  equa¬ 
tion.  However,  the  present  approximation  does  not  give  sufficiently  accurate 
results  for  the  strains  and  stresses  in  the  shell,  and  the  approximation 
using  five  constants  must  be  employed  when  these  quantities  are  required. 


-26- 


d)  Approximation  Using  Three  Constants  -  Membrane  Effects  Only. 
n  4  0 

As  in  part  (b)  of  this  section,  a  membrane  approximation  may  be  used 
to  obtain  the  frequencies  of  sufficiently  thin  shells  in  the  lower  modes 
of  n.  The  frequency  equation  and  the  mode  shapes  for  the  membrane  shell 
are  obtained  by  setting  the  value  of  the  constant  k  equal  to  zero  in 
Eq.  (l-k8)-(l-5l).  The  determinantal  frequency  equation  is  given  on  Page  2k. 

The  remarks  of  part  (c)  of  this  section  regarding  the  validity  and 
applicability  of  three  constant  approximation  also  hold  for  this  case<> 


e)  Mode  n  =  0 


It  is  shown  in  Reference  (l)  ^  that  the  exact  solution  of  the  bound¬ 

ary  value  problem  for  the  free  vibrations  of  a  thin  cylindrical  shell  in 
the  mode  n  =  0  yields  an  infinite  number  of  frequency  numbers  lying 
between  closely  spaced  limits.  For  a  steel  shell  (  \)  =  0.3)  with  fc  *  0.8, 
these  limits  are  given  by  H  =  0.8833  and  0.91*  Since  the  frequencies  and 
their  corresponding  sub-modes  are  so  close  together,  a  forced  vibration 
analysis  for  the  response  of  the  shell  in  the  mode  n  »  0  would  require  an 
expansion  involving,  in  general,  many  of  thpse  sub-modes  of  frequency  o>qj . 
This  is  not  practicable  and  in  general,  other  methods  must  be  employed  to 
approximate  the  response  of  the  shell  in  the  zero  mode.  This  problem  is 
discussed  in  some  detail  in  Section  VI  of  this  Report.  One  suitable 


method  is  to  consider  the  shell  as  a  series  of  separate  rings. 
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II  Comparison  of  Results  for  Membrane  Shells 

The  values  of  the  frequency  number  Mq,  and  the  normalized  strains  ae^,  aezz> 

obtained  from  the  five  constant  membrane  approximation  of  Section  I, 

Part  b,  will  be  compared  with  the  results  of  the  membrane  analysis  of 

ll  X 

Reference  (l).  An  empty  steel  shell  in  which  *=  0.8,  >*=  0.3  and  —  » 
will  be  considered.  For  the  modes  1  <  n  <  9>  the  shell  acts  as  a 
membrane,  as  shown  in  Fig.  (III-1). 

Table  (il-l)  gives  the  computed  values  of  the  lowest  frequency 
number,  M^,  for  the  modes  1  <  n  <  6  of  the  steel  tank  under  consideration* 

The  results  of  both  the  five  constant  and  the  three  constant  approximations 
are  given  and  these  results  are  compared  with  the  values  of  from  the 
membrane  theory  of  Reference  (l).  In  addition,  the  percentage  error  in 
the  frequency,  c 0^,  is  given  for  the  results  of  both  approximations. 

It  is  seen  that  the  error  in  the  frequency  varies  between  1.6$ 
and  6.3$  for  the  five  constant  approximation,  and  between  4.3$  and  19*1$ 
for  the  three  constant  approximation.  It  is  concluded  that  both  approxi¬ 
mate  theories  give  sufficiently  accurate  estimates  of  the  frequencies  of 
vibration  of  the  shell  for  use  in  problems  concerning  the  response  of  the 
shell  to  dynamic  loading.  It  should  be  noted  that  the  increased  accuracy 
of  the  five  constant  approximation  is  obtained  at  the  cost  of  the  extra 
labor  involved  in  expanding  a  fifth  order  determinant  rather  than  one  of 
the  third  order. 

For  applications  in  which  the  values  of  the  strains  and  stresses 
of  the  shell  in  each  mode  n  are  required,  the  five  constant  theory  must 
be  used.  It  should  be  noted  that  although  the  frequencies  of  vibration 


and  the  shape  of  the  modes,  u(0,  %),  v(0,  t)  and  w(0,  t)  can  be  estimated 
by  means  of  the  three  constant  theory,  the  strains  and  corresponding 
stresses  which  Involve  derivatives  of  the  displacements  cannot  be  ob¬ 
tained  from  this  theory  with  sufficient  accuracy  for  use  In  forced 
vibration  problems. 

On  the  other  hand,  the  five  constant  theory  gives  sufficiently 
accurate  values  of  the  shell  strains  and  stresses  in  the  range  of  parti* 
cuiar  Interest.  For  future  applications,  the  longitudinal  stress  o 
at  the  bottom  of  the  shell  will  be  of  particular  interest  in  analysing 
the  uplifting  of  the  shell  under  dynamic  loading.  Also  of  interest  will 
be  the  direct  stress  oQd  at  the  top  of  the  shell  for  possible  use  in 
connection  with  both  buckling  considerations.  The  determination  of  thesa 
stresses  requires  the  values  of  the  strainscgg  and  «  at  the  top  and 
bottom  of  the  tank. 

Fig.  (lI-l)-(lI-4)  show  the  variation  of  the  strain  ac  with  the 

%% 

height  of  the  tank  for  the  modes  n  ■  1,  2,  k  and  6.  The  strains  obtained 
from  both  the  five  constant  theory  and  the  membrane  analysis  of  Reference 
(1)  are  shown.  It  is  felt  that  the  results  of  the  five  constant  theory 
are  of  sufficient  accuracy  for  use  in  applications  to  forced  vibration 
problems. 

Fig.  (lX-5)-(II-Q)  show  the  variation  of  the  strain  a «0g  with  the  height 
of  the  tank  for  the  modes  n  ■  1,  2,  4,  6,  Again,  it  is  felt  that  the  results 
of  the  five  constant  theory  are  of  sufficient  accuracy  for  use  in  applica¬ 
tion  to  forced  vibration  problems. 

Fig.  (XI-9)- (11-12)  give  the  values  of  the  shear  strain  a«sg  for  the 
modes  n  ■  1,  2.  k  and  6.  It  is  seen  that  the  five  constant  approximation 


gives  a  poor  result  for  this  strain.  However,  as  this  strain  is  not  of 
any  practical  importance  in  our  applications,  this  is  no  deterrent  to  the 
use  of  the  five  constant  theory.  A  better  value  of  a  ez0  could  be  ob¬ 
tained  from  a  six  or  possibly  seven  constant  theory,  but  this  would  lead 
to  higher  order  frequency  determinants  and  would  serve  no  useful  purpose 
in  future  applications. 


Table  II-l 

Frequency  Humber, 

M 

n 

Empty  Steel  Tank  L/a 

*  0.8 

)) 

-  0.3 

h 

1 

a 

“  1200 

Five 

Constant 

Three 

Constant 

Reference  (l) 

Approximation 

Approximation 

n 

M  , 
nl 

Error  in 
%!<*> 

M 

nl 

Error  in 

V  W 

Mnl 

1 

.4575 

3.9 

.4610 

4.3 

.4242 

2 

.1906 

2.9 

.197^ 

4.7 

.1802 

3 

.0859 

1.5 

.0934 

5.8 

.0834 

4 

.0432 

1.6 

.0494 

8.8 

.0418 

5 

.0241 

3.4 

.0290 

13.3 

.0226 

6 

.0147 

6.3 

.0184 

19.1 

.0130 

0  M  .Eh 

2  nl 

111 '  VW) 

>S’VN\n**n»S':n“Cvv 

WV'-A^v 

s'  s’  •  ■  *.’ 

.  %  ■•*.»*«  %  •- 

.  V  *  •  ‘  *  -  -  *  k  * 

NORMALIZED  MEMBRANE  STRAINS -Q6 


NORMALIZED  MEMBRANE  STRAINS 


FIG. IE-3  NORMALIZED  MEMBRANE  STRAINS 


FiG.E-4  NORMALIZED  MEMBRANE  STRAINS-Q8 


FIG.IL-5  NORMALIZED  MEMBRANE  STRAINS 


NORMALIZED  MEMBRANE  STRAINS 


FlG.JC-7  NORMALIZED  MEMBRANE  STRAINS 


NORMALIZED  MEMBRANE  STRAINS 


FIG.ir-10  NORMALIZED  MEMBRANE  STRAINS 


NORMALIZE  MEMBRANE  STRAINS 


NORMALIZED  MEMBRANE  STRAINS 
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III  Bending  Effects  on  the  Frequencies  and  Modes  of  Free  Vibrations 
of  Thin  Cylindrical  Shells. 

This  section  presents  a  study  of  the  effect  of  bending  strains  on  the 
frequencies  and  modes  of  free  vibrations  of  thin  cylindrical  shells.  It 
should  be  noted  that  for  shells  which  are  sufficiently  thin  to  be  considered 
membrane  shells,  the  frequencies  of  free  vibrations  are  independent  of  the 
thickness  to  radius  ratio,  h/a,  and  depend  only  on  the  height  to  radius 
ratio,  (  *  L/a.  For  a  shell  of  a  given  material  and  ratio  £  lying  in 
this  range,  the  membrane  frequency  would  be  the  same  regardless  of  the 
■thickness  of  the  shell.  The  problem  therefore  takes  on  two  aspects  for 
discussion: 

1.  The  establishment  of  ranges  in.  which  shells  of  practical  interest 
could  be  considered  as  membranes. 

2.  The  derivation  of  a  simple  procedure  which  would  give  an  estimate 
of  the  effect  of  bending  strains  on  the  lowest  membrane  frequency  in  any 
mode  n.  Such  a  procedure  would  provide  a  Bimple  method  for  establishing 

the  range  of  modes  in  which  a  given  shell  could  be  considered  to  be  a  mem¬ 
brane.  It  would  also  provide  an  estimate  of  the  effect  of  bending  strains 
on  the  lowest  membrane  frequency  in  any  mode  n,  in  the  range  in  which  bend¬ 
ing  effects  are  of  importance.  The  advantage  of  obtaining  such  information 
without  the  necessity  of  expanding  the  determinantial  frequency  equation, 
aq.  (1-23),  for  each  shell  with  a  different  h/a  ratio  is  obvious. 

It  should  be  noted  however,  that  if  values  of  the  strains  and  stresses 
in  the  shell  are  also  required,  the  determinants!  equation,  Eq.  (1-23), 
must  be  expanded  to  obtain  the  lowest  root  M^,,  and  the  corresponding  mode 
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shspes  must  be  computed  from  Eq.  (l-l6)-(l-20).  The  use  of  the  mode 
shapes  resulting  from  membrane  theory  does  not  give  an  accurate  deter¬ 
mination  of  the  strains  in  the  shell  in  the  range  in  which  bending 
effects  are  of  importance. 

An  approximation  for  the  determination  of  the  frequencies  and  mode 
shapes  of  thin  cylindrical  shells  with  small  —  ratios  was  given  in 

£t 

Section  I,  parts  (b)  and  (d),  by  considering  the  shell  to  be  a  membrane 
with  no  bending  stiffness.  The  results  of  such  an  approximation  were 
noted  to  be  valid  for  modes  with  a  low  circumferential  wave  number  n. 

In  order  to  establish  a  range  of  validity  of  -  and  n  for  this  approxi- 

£1 

nation  ^  and  to  study  the  affects  of  bending  strains  on  the  frequencies 
of  the  shell,  computations  were  made  and  are  presented  for  the  following 
two  shells  of  practical  interest! 

a)  Unprotected  steel  tank,  -  »  “  ■  0.8,  V  *»  0.3 

b)  Steel  tank  with  concrete  shielding,  -  ■  -p-  ,  ^  »  0.8,  ^  -  0. 3 

Figure  (III-1)  shows  a  curve  of  the  frequency  number  Mn  (correspond¬ 
ing  to  the  lowest  frequency  j  a  1  in  each  mode  n)  plotted  against  the  mode 
number  n  for  the  unprotected  steel  tank.  The  values  of  on  the  curve 
labeled  "membrane  only"  were  derived  from  the  membrane  frequency  equation, 
Iq.  (1-41),  while  those  on  the  curve  labeled  "membrane  plus  bending"  were 
derived  from  the  complete  frequency  equation,  Eq.  (1-23).  For  a  shell 
of  this  relatively  small  thickness  to  radius  ratio,  both  theories  give. the 


(1)  This  range  is  quite  different  for  Bteel  tanks  and  for  steel  tanks 
with  concrete  shielding. 


same  frequencies  in  the  ranee  1  <  n  <  9*  For  the  modes  in  which 
9  <  n  <  14,  the  inexbensional  bendins  strains  and  the  extenoional  membrane 
strains  have  about  equal  effects  on  the  frequency.  For  n  >  14,  the 
bending  strains  become  predominant  and  control  the  frequencies  of 
vibration  of  the  shell. 

It  should  be  noted  therefore,  that  if  these  modes  are  to  be  used  in 
the  analysis  of  the  forced  vibrations  of  the  empty  steel  tank  to  dynamic 
loading,  the  membrane  theory  can  be  used  and  bending  effects  may  be 
neglected.  In  such  an  analysis,  practiced  considerations  of  the  convergence 
of  the  mode  series  for  the  response  of  the  shell  to  dynamic  loadings,  would 
certainly  not  require  more  than  the  number  of  accurate  modes  given  by  the 
membrane  analysis. 

The  situation  changes  radically  when  the  steel  tank  1b  protected  by 
a  concrete  shielding.  Fig.  (III-2)  shows  the  curve  of  Mn  versus  n  for 
this  case.  As  in  Fig.  (Ill-l),  the  "membrane  only"  and  "membrane  plus 
bending"  curves  are  computed  from  Eq.  (l-4l)  and  Eq.  (1-23)  respectively. 
Since  the  thickness  to  radius  ratio  of  the  protected  shell  is  roughly 
forty  times  larger  than  that  for  the  unprotected  shell,  it  is  to  be 
ejected  that  the  bending  effects  will  become  of  importance  at  a  much 
lower  mode  n  than  in  the  case  of  the  unprotected  shell.  It  is  seen  from 
Fig.  (III-2)  that  a  membrane  analysis  gives  accurate  results  in  the 
range  1  <  n  <  3*  Bending  effects  start  to  become  of  importance  at  n  *  3 
and  for  n  >  4,  the  bending  strains  become  predominant  and  control  the 
frequency  of  vibration  of  the  shell.  The  forced  vibration  analysis  of 
the  response  for  a  protected  steel  tank  to  dynamic  loading  would  therefore 
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re  quire  that  bending  effects  be  included  in  all  modes  where  n  >  3* 

A  simple  estimate  of  the  effect  of  bending  strains  on  the  membrane 
frequencies  of  a  cylindrical  shell  will  now  be  developed.  This  correction 
to  the  membrane  frequency  will  enable  the  computation  of  the  frequencies 
of  vibration  of  shells  in  the  range  where  bending  energy  is  of  importance, 
without  resorting  to  an  expansion  of  the  fifth  order  determinantal  fre¬ 
quency  equation,  Eq.  (1-23),  for  each  shell  with  a  different  h/a  ratio. 
'Moreover,  for  a  shell  with  a  particular  thickness  to  radius  ratio  h/a, 
it  will  enable  us  to  determine  the  range  in  which  the  shell  may  be  con¬ 
sidered  to  be  a  membrane. 

To  obtain  this  bending  correction  to  the  membrane  frequencies  of  a 
shell,  consider  the  case  of  purely  inextensional  motions  of  a  thin 
cylindrical  shell  in  free  vibrations.  The  condition  of  inextension  of 
the  middle  surface  of  the  thin  shell  requires  that  the  strain  ac^  *» 

Vg  -  W  be  equal  to  zero.  This  condition  is  satisfied  by  displacements 
of  the  form 

v(z,  0,  t)  =»  |  cosnS  e1^  (ill-l) 

v(z,  t)  »  £  |  sinnS  e1^1  (III-2) 

u(z,  0,  t)  =  -|  cosn0  e1<uTt  (III-3) 

n 

The  frequency  of  these  inextensional  oscillations  can  be  determined 
by  Rayleigh's  Principle: 


T 

max 


V 

max 


(m-4) 


Substituting  Eq.  (III-l)-(lII-3)  into  the  strain  energy  expression  for 
V£  as  given  by  Eq.  (1-12)  the  potential  energy  of  the  shell  is  given  by 


2raax 


Eh* 

24(1-V2) 


(XXX. 5) 


while  the  maximum  kinetic  energy  le  obtained  by  substitution  of  these 
equations  into  Eq.  (1-7) : 


T 

max 


2  2 
mina  tttj 


(in-6) 


Substituting  Eq.  (III-5)  and  (III-6)  into  Eq.  (III-4)  the  frequency  number 
ia  given  by  the  relation 


..  h2  n2(l-n2)2 

Mnl  "  7T5  2“ 

12a  n  ♦  1 


2.2 


2(1-1?)  +jQ‘ 

JL 


2.2 


n2+l 


(XXX-7) 


where  i*  defined  by  Xq.  (1-24) .  The  corresponding  frequency  ^  becomes: 


2 

%1 


Eh 


12 


fid-if)* 


n2(ng-l)2 

n2+l 


1+ 


JZkgl 


1  + 


<n2»l)n2ia 


(Ill-ti) 


Xq.  (III-7)  and  (III-8)  give  the  frequencies  of  the  purely  inextensional 
motions  of  thin  cylindrical  shells. 


(2)  The  potential  energy  of  extension,  is  identically  equal  to  zero 


for  these  motions. 


To  obtain  an  estimate  of  the  effect  of  bending  strains  on  the  low¬ 
est  frequency  factor  in  any  node  n^  the  value  of  for  the  membrane 
shell  only  is  computed  from  Eq.  (l-4l)  or  Eq.  (1-57)*  Let  the  quantity 
be  called  M^.  The  value  of  the  frequency  factor  for  purely  in- 
extenoional  motions  of  the  shell  is  computed  from  Eq.  (III-7).  Noting 
that  the  factors  M  are  proportional  to  the  square  of  the  frequency,  an 
approximate  value  of  the  frequency  factor  which  includes  both  membrane 
and  bending  effects  is  obtained  from  the  relation: 


M  =  M  +  M  _ 
n  nra  nl 


(m-9) 


The  corrected  frequency  is  then  given  by  the  relation 


f  2  2 

00  =  0)  +  (0  T 
n  1  am  nl 


(m-10) 


whore  (oirj  and  a>nI  are  confuted  from  Eq.  (1-24)  and  Eq.  (III-8)  res¬ 
pectively.  The  method  used  in  obtaining  Eq.  (III-10)  iB  sometimes  known 
as  Southwell's  method. 

Tables  (ill-l)  and  (III-2)  show  the  application  of  Eq.  (III-9)  to 
the  cases  of  the  protected  and  unprotected  steel  tanks  respectively. 
Column  4  of  each  table  gives  the  value  of  Mq  computed  from  Eq.  (IIX-9) 
while  column  5  gives  the  value  obtained  by  expanding  the  fifth  order 
determinontal  frequency  equation,  Eq.  (1-23).  It  may  be  noted  that  very- 
good  approximations  to  the  frequency  factor  can  be  obtained  using 
Eq.  (III-9).  In  addition,  the  range  of  applicability  of  the  membrane 
theory  can  easily  be  established  from  such  tables,  thus  verifying  the 
results  shown  in  Fig.  (Ill-l)  and  Fig.  (III-2). 

It  should  again  be  noted  however,  that  when  strains  and  stresses  are 


required  in  a  mode  in  which  bending  effects  are  important,  the  value  of 
should  be  obtained  by  expanding  Eq.  (1-23)  and  the  correct  mode  shapes  should 
be  evaluated  using  Eq,  (l-l6)-(l-20). 


Table  III-S 


Effect  of  Bending  on  Membrane  Frequencies 


Empty  Steel  Tank  -  No  Shielding 
Five  Constant  Approximation 
L/a  -  0.8,  hi  V)  -  0*3 
1  '  ~  *  Tpoo* 


a 


Mnl 

% 

Iq.  (III-9) 

n 

Kq.  (1-23) 

.00004 

.02416 

.02412 

.0001 

.0148 

.0147 

.00060 

.00420 

.00420 

.00123 

.00345 

.00350 

.00298 

.00422 

.00421 

.0093^ 

.00995 

.00990 

IV  Free  Vibrations  of  Fluid  Filled  Shells 


Determination  of  the  Virtual  Mass  of  the  Fluid 


Methods  leading  to  the  determination  of  the  frequencies  and  inodes  of 


free  vibration  of  empty  cylindrical  fuel  tanks  have  been  presented  in  pre* 


vious  sections  of  this  report.  In  this  section,  the  analysis  will  be 


extended  to  the  case  of  partially  full  cylindrical  fuel  tanks.  As  in  the 


case  of  the  empty  tank,  the  modes  of  vibration  will  be  characterized  by 


the  parameter  n,  the  number  of  circumferential  waves  in  the  mode. 


Consider  the  shell  of  Figure (tV-l)  which  is  filled  to  a  height  yl 


with  an  incompressible  inviscid  fluid  of  density  p  ..  The  assumption  of 


the  incompressibility  and  zero  viscosity  of  the  fluid  implies  the  existence 


of  a  velocity  potential  function,  0  (r,0,  z,  t),  such  that 


*  ?£  ♦  4  rl  -  r§  -  o 

dr2  ^  r2  d02  dz2 


(IV-1) 


The  motion  of  the  fluid  may  be  derived  from  0,  since  the  velocity  of  the 


fluid  in  a  direction  s  is  given  by  the  directional  derivative 


(IV-2) 


The  expressions  for  the  radial,  tangential  and  longitudinal  velocities  of 


the  fluid  particles  in  polar  coordinates  are 


vr  "  *  5? 


(IV-3) 


>8  -  *  S 


(IV-4) 


-s 


(IV-5) 
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The  pressure  in  the  fluid  at  any  point  in  space  is  given  by  the  deriva¬ 
tive  of  0  with  respect  to  time, 


P  - 


(IV-6) 


where  it  has  been  assumed  that  the  velocities  are  small. 

The  motion  of  the  fluid  during  the  free  vibrations  of  the  partially 
filled  shell  can  be  obtained  by  the  solution  of  the  boundary  value  problem 
involving  Eq.  (iV-l)  and  the  following  three  boundary  conditions: 

a)  The  radial  velocity  vr  of  the  fluid  must  be  equal  to  the  radial 
velocity  w  of  the  shell  on  the  surface  r  -  a. 


w(t,  6,  t)  - 


or 


Jr  »  a 


(IY-7) 


b)  The  longitudinal  velocity  v  at  the  bottom  of  the  shell,  z 
must  be  equal  to  zero. 


z  oz 


z  « 


.  Ik 


(IV-8) 


and 


c)  On  the  free  surface  at  x  ■  +  ,  the  pressure  must  be  equal  to 


zero: 


_  80(r,  Q,  z,  t) 

p  .  P  '  -t 


0  . 


7k 


(lv-9) 


Eq.  (IV-l),  (IV-7),  (IV-8),  (IV-9)  enable  the  determination  of  the  velocity 
P'  ’jnttal  0  and  subsequently  of  the  kinetic  energy  of  the  fluid  from  the 


relation: 


fluid 


(IV-10) 


where  the  integration  la  taken  over  all  surface*  of  the  fluid  and  is 
the  normal  derivative  of  0  on  the  particular  surface  over  which  the 
Integration  is  performed.  The  kinetic  energy  may  also  be  obtained  by  a 
volume  integration  of  the  square  of  the  velocity* 


Tf  -  |  JJJ  til  *  4  +  rdrd6dt 


(lV-11-a) 


which,  upon  substitution  of  Eq.  (IV-3)  •  (IV-5 )#  becomes 


*  ■  *  & 


2  +  <M>2  +  4  ®2  rdrd0d* 

r 


(lV-11-b) 


The  determination  of  0  will  take  place  in  three  steps  using  the 
principle  of  superposition.  Three  separate  potentials,  0^,  0g,  and  0^ 
will  be  obtained  such  that 


0  -  +  02  +  ^3 


(IV-12) 


will  satisfy  the  boundary  value  problem. 

The  various  potential  functions  0^(r,  6,  z,  t)  of  Iq.  (IV-12)  can  be 
expressed  as  a  summation  of  their  respective  components  in  the  n  modest 


0k  "  ^  6 ’  *' 


(IV-13) 


[1)  tf'or  example,  "Hydrodynamics 
L932,  pp.46,  Eq.  (4). 


"  by  H.  Lamb,  Dover  Publications,  Sixth  Edition, 
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In  the  analysis  that  follows,  the  expression#  for  will  he 

determined  and  th#  virtual  mass  m  in  each  mode  will  be  evaluated. 

vn 

Expressions  for^  ^ ,  the  displacement  of  the  fluid  surface  will 
also  be  presented.  It  i#  convenient  to  consider  the  problem  for  the 
following  three  cases:  a)  n  0  0  -  Three  Constant  Approximation, 
b)  n  «  0  and  c)  n  /  0  -  Five  Constant  Approximation. 


(a)  n  4  Q  .  Three  Constant  Approximation 


Consider  the  potential  function 
r  In 

*tn  ■  ’  Bn  I  +  d  “S3  C0‘  nS  (IV'1',) 

J  na 

p 

which  satisifes  the  equation^  0-^  *  0.  The  component  of  velocity  in 
the  radial  direction  is 

vrl  ■  Bn[l  +  d]  <S)"’X  «•  »•  «.<*)  *  <1V-l5> 

Equating  the  radial  velocity  of  the  shell  to  the  radial  velocity  of 
the  fluid  at  r  ■  a,  Eq.  (IV-7)  nay  he  used  to  evaluate  Bq.  Defining  the 
generalized  coordinate  of  the  n  shell  mode  as  q^t)*  the  shell  displace* 
meats  are  given  by  the  relations 


u  -  u(z,  0)qn(t) 

(1V-16)  (a) 

v  -  v(z,  0)qn(t) 

(b) 

w  -  w(z,  0)qn(t) 

(c) 

where  u(s,  0),  v(a,  0)  said  w(s,  0)  are  defined  by  Eq.  (I-53)-(I«55)»  The 
value  of  v  is  obtained  by  differentiating  Eq.  (lV-l6-c)  and  shifting  the 
axis  of  s  coordinates  to  that  of  Fig.  (IV-1): 

w  -  -  °n  *£  [§  *  d]  00‘  aS  V(t)  *  (1T’1T) 
Substituting  Eq.  (IV-15)  and  (IV-17)  into  (IV-7),  the  value  of  Bn  Is 
Bfl  »  -  2jjjf  ,  and  the  velocity  components  due  to  0^  becomes 

vrl  ’  •  °n  ^  [l  +  d\  (a)D'1  co*  nS  V(t)  (IT-l8) 


(2)  The  minus  sign  in  Eq.  (IV-17)  is  required  to  refer  the  shell  velocity 
to  a  "positive  velocity  in  the  positive  r  direction"  convention  used  in 
Eq.  (IV-2)  and  (IV-7). 


The  fluid  potential  0ln  gives  a  longitudinal  velocity  component  v^ 
which  is  independent  of  z.  In  order  to  satisfy  the  boundary  condition 


of  Eq,  (IV-8)  a  second  function  0gn  muBt  be  added  to  0^. 

p 

The  velocity  potential  0gn(r,  Q,  z,  t)  is  chosen  so  that  <?0gQ  »  0 
and  that  in  conjunction  with  0^,  the  boundary  condition  of  Eq.  (IV-8)  is 
satisfied.  At  the  same  time,  the  radial  velocity  Vrg  is  chosen  so  that 


r*a  r«a 


(IV-21) 


since  already  satisfies  the  boundary  condition  of  Kq.  (IV-7). 
The  form  of  0g  is  taken  as 

02n  -  Y  °ni  8iDh  ^ST)  Jn(“T“)  cos  (IV-22) 

1»1 


where,  in  order  to  satisfy  Eq.  (IV-21),  the  aQ^  are  the  zeros  of  the  deri 
vative  of  the  Bessel  Function  of  order  n,  1a 


J»  (a  4 ) 
n  'ni' 


(IV-23) 


The  component  of  velocity  in  the  longitudinal  direction  is 


0-3 


C-0$h  (*n‘£) 


CL 


) COS  ft# 


C=t 


/f  , 

(7Z-24) 


The  boundary  condition  of  Eq.  (IV-8)  may  now  be  written  as 


[V ° 

z 

and  upon  substitution  of  Eq.  (IV-202  and  (IV-24),  Eq.  (IV-25)  becomes 

Cosh^L)J-n(«^r)  =^n^nan-,  0*-*fc) 

The  constant  C  .  can  be  evaluated  using  the  orthogonality  property 
nl 

of  Bessel  Functions.  Multiplying  both  sides  of  Eq.  (1V-26)  by  the 
a  ,r 

function  where  is  any  one  of  the  roots  of  (a^)  *  0, 

end  integrating  with  respect  to  r,  one  obtains 


The  orthogonality  relation  for  the  Bessel  functions 


(S) 


is 


(3)  Bee,  for  example,  McLachlan,  N.W.,  "Bessel  Functions  for  Engineers", 
Oxford  University  Press,  1934,  Pg.  160,  Eq«  (47)  and  (48). 


Also  used  will  be  the  relation 

n  it • 

which  is  valid  for  J*(a  .)  »  0. 

n'  ni7 

Using  Eq.  (IV-28)  and  (IV-29),  the  constant  is  evaluated  from 
Eq.  IV-2T): 


(iz-ze) 


The  velocity  components  due  to  the  function  become: 


^2  Z  cm  ^/nh(^f^)Jn  cos  ns  Cfn (t) 

vez~  Z  £/>l  ■s'nh(^£)^n sm  ne p(*) (jr-zz) 


V«2  which  is  given  by  Eq#  (IV-24) 


The  potential  function  obtained  by  adding  0ly|  +  0gQ  satisfies  (Eq.  IV-l) 
and  represents  a  motion  in  which  the  boundary  conditions  of  Eq.  (IV-7)  and 
Eq.  (IV-8)  are  satisfied.  Moreover,  this  motion  maintains  a  plane  surface 
at  the  top  of  the  fluid,  z  ■  on  which  the  particles  have  zero  velocity. 

To  satisfy  the  condition  of  zero  pressure  on  the  free  surface,  a  third 
fluid  potential  0^n(r,  6  >  z>  ^  *8  introduced.  This  function  represents 
the  sloshing  motions  of  the  fluid  on  a  partially  filled  tank  under  the 
conditions  of  zero  velocity  at  the  tank  bottom  and  zero  radial  velocity 
at  the  surface  of  the  tank  r  =  a.  Thus,  it  represents  a  motion  in  which 
the  tank  remains  rigid  while  the  fluid  inside  it  moves.  The  velocity 
potential  0^n (ri  G  t)  is  chosen  as 


where  the  values  of  are  again  the  zeros  of  the  derivative  of  the 
Bessel  Function  of  order  n  as  given  by  Eq.  (IV-23).  The  coordinates 
A^(t)  could  be  considered  as  the  generalized  coordinates  of  the  sloshing 
modes  of  the  fluid  which  are  to  be  superimposed  on  the  motions  due  to 
the  tank  displacements  given  by  the  potential  functions  +  0^.  The 
velocity  components  due  to  the  function  0^  are 

Vr3  =  1  Am  cos  „eti(°^)coshl f-jj 

W  Cm  34) 


A 


I 
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A nL  na  s/n  Jn  C  osh\j^r(z  -r^j/ 

hi  r<*nL  Jn  (o<ni)*tnh  (jJT-3S) 


rt 


3 


3 


a 

s 

a 


and 


AW-  cos  tH 

J^fanu)  smh  (^r)L  (pr-  3Cf) 

Eq*  (lV-3^)  and  ( 17-36)  show  that  *  0  and  that  ^  »  0. 

The  coordinates  AQi(t)  are  as  yet  undetermined  but  may  be  expressed  in 
terms  of  the  ^(t),  the  generalised  coordinates  of  the  shell  motion.  This 
is  done  by  utilisation  of  the  condition  that  the  fluid  surface  at  s  • 
must  be  a  free  surface  with  tero  pressure. 

To  find  this  relation,  consider  the  general  fluid  motion  governed 
by  the  potential  function  0  of  Bq.  (IY-13)  so  that  the  velocity  of  the 
fluid  in  any  direction  •  is  given  by 


V2" 


rfll 


+  Vs2  +  Vs3 


(IV-37) 


Hie  kinetic  energy  of  the  fluid  may  be  determined  from  Iq.  (lV-10), 
which  upon  substitution  of  the  appropriate  values  of  and  its  space 
derivatives  on  the  surfaces  r  «  a  and  s  ■  +  becomes j 
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Kubotituting  the  values  of  0n  end  its  ppace  derivatives  Into  Kq.  (IV- 
38),  and  performing  the  required  integrations,  the  kinetic  energy  of  the 
fluid  on  the  nth  mode  becomes: 


Eq.  (IV-39)  can  be  written  ao 


-88- 


vhere 


i»v 

Let  the  kinetic  energy  of  the  shell  in  the  n  code  be  given  as 


^nshell 


(IV- 45) 


and  the  potential  energy  of  the  shell  he 


^nahell 


(IV.46) 


The  equations  of  motion  governing  the  generalized  coordinates  q^(t)  and 
A^(t)  may  he  obtained  wy  the  use  of  Lagrange s  equations: 


-  (H ) 

dt 


&v 

4*  r  ■■■ 
OQ, 


V  dqk 


(rv-47) 
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vhere  T  la  the  total  kinetic  energy  of  the  system  in  the  x  “  mode  and  is 
given  by  the  sum  of  £q.  (IV-40)  and  (IV-45);  V  is  the  potential  energy  of 
the  system  and  is  ^iven  by  Eq.  ( IV— 46) j  and  is  the  generalised  force  in 
the  mode  k  due  to  any  external  loading  and  to  body  forces  such  as  the 
force  of  gravity. 

K  L 

Noting  that  the  pressure  on  the  fluid  surface  t  »  is  tero  and 
neglecting  the  effect  of  the  gravity  forces^,  the  generalized  force 
is  zero  for  free  vibrations  and  the  Lagrange  equations  on  the 
generalized  coordinates  and  become: 


(2-D+Mn)j„  +(2lpd)%»  +  +Ki$"'a0 


•  0 


uzi  p  +  -G 


/ . 


(  t=  1,1,3  ■■  '  ~j 


(jS-4s) 


vhere  Eq.  (IV-49)  holds  for  each  value  of  1. 


(4)  If  gravity  is  included  in  the  analysis,  there  is  a  contribution  to 
the  generalized  force  term,  Q^,  vhich  contains  terms  A^.  The  analysis 
of  this  complete  system  gives  rise  to  two  sets  of  modes,  one  of  vhich  is 
very  close  to  the  tank  modes  qn  which  ere  obtained  when  gravity  is  neg¬ 
lected.  The  second  set  of  modes  are  very  nearly  the  sloshing  modes  of 
the  fluid  in  a  rigid  tank  in  the  presence  of  gravity.  These  modes  have 
extremely  long  periods  as  compared  to  the  qn  modes  and  give  practically 
no  contribution  to  the  displacements  and  stresses  in  the  tank.  The  grav¬ 
ity  sloshing  modes  are  required  for  the  analysis  of  possible  fluid  spill¬ 
ing  from  the  tank  and  will  be  considered  in  a  later  section  of  the  report. 

The  modes  qfl  obtained  by  neglecting  gravity  will  thus  be  used  to  determine 
the  stresses  and  displacements  in  the  shell. 


Solving  Eq.  (IV-49)  for  Afli, 

Ar?l  “  "  Bm  (i*  °*)  Q&-So) 

20 si  * 

and  substituting  Eq.  (IV-50)  into  Eq.  (iV-to),  we  obtain  the  equation 
on  the  generalized  coordinate  q^t), 


Equation  (IV-51)  implies  that  the  kinetic  energy  of  the  fluid,  Eq.  (IV-39), 
(IV- to)  can  then  he  written  in  terms  of  the  generalized  coordinate  q^t) 
only: 

^//w  **  £'>$»(*)  (jz-sz) 

where  the  value  of  CQ  is  given  by 

^cx-i^-ekj  firs,) 


Substituting  the  approximate  constants  into  Eq.  (IV-52)  and  Eq.  (IV-53), 
the  kinetic  energy  of  the  fluid  in  the  nth  mode  becomes: 


The  "virtual  mass"  of  the  fluid  in  the  n  mode,  m  .  is  defined  so 

vn 

iy. 

that  the  kinetic  energy  of  the  fluid  in  the  n  mode  is  given  by 


,Z'!i Z2L 

=  /  /  wffe  0,  t)  aSSo/zt 


(zr-sej 


vhere  vn(z ,  &  ,  t),  the  radial  displacement  of  the  shell  is  given  by 
Eq.  (IV-l6c)e  The  total  kinetic  energy  of  the  partially  filled  shell 
plus  the  fluid  in  the  n °  mode  becomes 


(7*fL  ,2V f*L 

Jo  Jo  '• 


The  virtual  mass  of  the  fluid  can  thus  be  considered  to  be  an  additional 
tank  mass  moving  only  in  the  radial  direction. 

Substituting  Eq.  (lV-l6c)  into  Eq.  (IV- 55)  and  equating  the  result 
to  Eq.  (IV-54),  the  virtual  mass  of  the  fluid  in  the  mode,  n^, 
becomes: 


Lit* fit. 


*  *■ 


ZU»,  irL 


For  convenience  Eq.  (IV-57)  way  be  written  as 


m. 


vn 


V* 


(IV-5B) 


where  6^  is  the  quantity  inside  the  braces.  Figure  (IV-2)  shows  the 
variation  of  en  with  the  mode  n  for  the  various  fillings  denoted  by 
7  *  0.5,  0.75,  0.9,  and  1.0.  The  computations  were  made  for  a  tank  for 
which  L/a  «  0.8  and  include  terms  up  to  and  including  i  «  3  in  the  summa¬ 
tion  in  €  . 

n 

The  motion  of  the  fluid  surface  at  z  ■  +  in  each  mode  n,  is 
required  for  the  analysis  of  the  response  of  the  partially  filled  tank  to 
dynamic  loading.  This  analysis  is  given  in  Chapter  VI  of  this  Report. 
These  motions  can  be  obtained  by  an  integration  with  respect  to  time  of 
the  longitudinal  velocity  V  on  the  fluid  surface  z  m  +  Zk  ,  Using 
Sq.  (IV-37)  and  noting  that 

(IV-59) 


the  velocity  of  the  fluid  on  the  surface  z  -  +  is  given  by  V  _.  Usinf 

^  Zj 


Eq.  (IV-36)  and  Iq.  (IV-50),  the  velocity  becomes: 


1  ]  =  SlCn  ill 

Z*n<.!jL\fonh(i 


(zr- 


Integrating  with  respect  to  time  and  noting  that  the  initial  displacement 
of  the  fluid  is  zero,  the  fluid  displacement  in  the  n^  mode  (Measured 
positive  upward), 1  (r,  Q,  t)  ie  given  by  the  relation 


-93- 


where 


Rewriting  Eq.  (lV-6lb)  as 


*  ni^1  6>  “  *T“  PniJn  008  n®  qn(t)  (lV-62a) 

where 

_  Z  °W  7 ‘arthfio^  §9  r  _/ _  /  / _ ^_/ _ 

P^-  "  X,fa»X*™-'/,V  U  (**.%) 


and  considering  a  tank  with  L/a  *  0.8,  the  values  of  the  coefficients  0^ 
are  given  in  Figures  (lV-3)-{lV-8)  for  the  modes  n  »  1-6  and  i  ■  1,  2,  3* 
In  each  case,  the  coefficient  is  given  for  tank  fillings  ranging  from 
y  -  0.50  to  7  -  1.00.  Table  (IV-1)  gives  the  corresponding  values  of 
the  roots  a  .  for  use  in  Eq.  (lV-62a). 


Consider  the  potential  function  0^  (r,  z,  t) 


0i"l  [r  m  (z  +  zi)2 

(IV-63) 

g 

which  satisfies  the  equation^/  0-^*0. 

The  component  of  velocity  in 

the  radial  direction  is 

vrl  “  “  a  * 

(IV-64) 

For  application  in  the  forced  vibration  analyai*  of  Section  (VI), a  t- 
indcpendent  radial  displacement  of  the  shell,  wq  ■  q^t),  vi31  be  con¬ 
sidered,  where  q(t)  ia  the  generalized  coordinate  of  the  shell  mode 
n  *  0. 


Noting  that  the  radial  velocity  of  the  shell  at  r 
.  (5) 
yo  “  *  % 


a  is  given  by 
(IV-65) 


the  boundary  condition  of  Eq.  (IV-7)  is  satisfied  by  the  potential 
function  0^  (r,  s,  t). 

The  velocity  conponenta  in  the  tangential  and  longitudinal  directions 


are  respectively 
ve  -  0 


(IV-66) 


\  "  +  I  *  2i)  5o<t)  '  (XV-6T) 

The  fluid  potential  0^  gives  a  longitudinal  velocity  component  vjE  «  0  at 
the  bottom  of  the  tank,  z  »  -  ,  thus  satisfying  the  boundary  condition 


(5)  See  footnote  (2)  on  Page  (IV-5) 


-95- 


of  T.q.  (IV-8). 

To  satisfy  the  condition  of  zero  pressure  on  the  free  surface,  a 
second  fluid  potential  02(r,  z,  t)  is  introduced.  This  function  repre¬ 
sents  the  sloshing  motion  of  the  fluid  in  a  partially  filled  rigid  tank. 
The  velocity  potential  0^  is  chosen  as 

6  (0i,ir) = -f  a*  *■& 

1  fa  d'o  (°<oL)$/f)h(o<otZ£) 


where  the  values  of  a  .  are  the  zeros  of  the  derivative  of  the  Bessel 

nl 

Functions  of  order  zero, 

Jo  <%i>  "  -  Ji  (aoi)  “  0  (IV-69) 

The  velocity  components  due  to  the  function  02  are 

KZ  JoC'Xoi)  (ZZ-7o) 

K52=  o 


(it- 72) 


Eq.  (IV-70)  and  (IV-72)  show  that  V 


- 

2 


■0  and  that  v 


21  r»a 


-*§ 


•  0.  The 


corodinates  AQ^(t)  may  be  expressed  in  terms  of  the  generalized  coordinate 
of  the  shell,  qQ(t)  in  a  manner  similar  to  the  case  where  n-  /  0. 


To  find  this  relation,  consider  the  general  fluid  motion  governed  by 
the  potential  function  0  =  0-,_  +  02  so  that  the  velocity  of  the  fluid  in 
any  direction  s  is  given  by  Vg  »  Vgl  +  Vs2»  The  kinetic  energy  of  the 
fluid  may  be  determined  from  Eq.  (IV-11), 


-t 


hj  (fii+  <Pzj]  I  rdrd&  dz 

■  '  *  (jZ-73; 


which  upon  substitution  of  the  appropriate  values  of  the  space  derivatives 
of  0  becomes: 

(lS-7*) 


Eq.  (IV-74)  may  be  written  as 

7T-W  =  ^OZ  +  Y  Do*  A  ol 

where 


and 


nnaJ  _ 

a  —C — - a  .  y\x 

2a  ,tunh(~ — ) 


(jr~7£r) 


(iv-76) 


(iv-77) 


\lv-78) 


Tofluld  -  Co  %  <*> 


(IV-05) 


where  the  value  of 


CO 


lal 


CQ  is  given  by 


(IV-86) 


fuhstituting  the  appropriate  constants  into  Eq.  (IV-85)  and  Eq.  ( IV-86), 
the  kinetic  energy  of  the  fluid  in  the  0  mode  becomes 


_  zS. 

z  .  V  Zo. 

L.c<c? 


The  "virtual  mass"  of  the  fluid  in  the  zero  mode,  mvQ,  is  defined 
so  that  the  kinetic  energy  of  the  fluid  is  given  by 


/ZW  (if L 

x„,.r  f3/  &-.>) 

and  wq  is  given  by  Eq.  (IV-65).  Substituting  this  value  of  wq  into 
Eq»  (IV-88)  and  equating  Eq.  (IV-87)  and  Eq.  (IV-88),  the  virtual  mass 
of  the  fluid,  becomes 


(nr- S3) 


The  cotton  of  the  fluid  surface  at  a  C  +  ~  In  the  zero  mode  Is 
obtained  by  an  integration  with  respect  to  time  of  the  longitudinal 
velocity  V  of  the  fluid  surface  z  »  +  ~  .  This  velocity  is  given  by 
the  sum  of  Eq.  (IV-67)  and  Eq.  (IV-72).  Using  the  relation  of  Eq.  (IV-83) 


the  velocity  becomes 


(zZ-<?o) 


Integrating  with  respect  to  time  and  noting  that  the  initial  displace¬ 
ment  of  the  fluid  is  zero,  the  fluid  displacement  at  the  surface  *  ■  +  -g 
in  the  mode  n  ■  0  (measured  positive  upward),"^  Q(r,  t),  is  given  by  ths 
relation 


»o 


'  iT/ 

(jT-Bt) 

uifiere. 

-1  (r,t)=  2  fB(t) 

(tt-BZ) 

Writing  Eq.  (IV-92)  as 

Vr*  *>  -  »oi  $ 

(IV-93) 

the  values  of  0Qi  are  given  in  Fig.  (IV-9)  for  a  tank  in  which  *  ■  0.8. 
The  coefficients  0Qi  are  given  for  tank  fillings  ranging  from  7  »  0.J0  to 
1,00  and  for  i  ■  1,  2  and  3* 

Using  Eq.  (IV-58)  with  n  «  0  and  (IV-89),  the  values  of  «Q  for  a 
tank  in  which  L/a  ■  0.8  are  shown  in  Fig.  (IV-2)  for  the  fillings  7  »  0.5, 


c)  a  /  0  Five  Constant  Approximation 


The  procedure  used  in  the  present  section  for  the  determination  of 
the  virtual  mass  and  the  fluid  displacements  is  exactly  analagous  to 
the  procedure  used  for  the  three  constant  approximation  of  Section  (IV-a). 


Consider  the  potential  function  0^Q 
**Ln  *  *  Bn  j°n  +  5Ej  “a-1  c0“e 


(IV-94) 


°n  - 1?  • « ♦  £  r 

n 


(IV-95) 


P 

which  satisfies  the  equation  y  0^  «  0 .  The  component  of  velocity  in 


the  radial  direction  is 


'rl  ‘  Bn  +  tTI  (f)”*1  '  (1V-96) 


Equating  the  radial  velocity  of  the  shell  to  the  radial  velocity  of  the 
fluid  at  r  ■  a,  Eq,  (IV-7)  may  be  used  to  evaluate  Bn»  Defining  the 

It 

generalised  coordinate  of  the  n  shell  mode  as  q^(t),  the  shell  displace¬ 
ments  are  given  by  Eq.  (IV-16)  where  u(t,  9),  v(s,6)  and  w(*,0)  are 
defined  by  Eq.  (l-25)-(l-27).  The  value  of  w  is  obtained  by  differentiat¬ 
ing  Eq,  (l-2?)  end  shifting  the  axis  of  z  coordinates  to  that  of  Fig.  (IV-1) 

*  "  *  cn  ^  L°n  +  rl]  COfln0  *  (*V-97) 

Substituting  Eq.  (IV-96)  and  (IV-97)  into  (IV-7),  the  value  of  B  is 

VL 


B  »  - 


(IV-98) 
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and  the  velocity  components  due  to  become 
Trl  •  *  Cn  ^  [  °n  +  ]  'I'0'1  coon9 

Vi"  (|>n"1  8lra9 


(IV-99) 

(IV-100) 

{TV-101) 


Proceding  as  in  Section  (V-a)  and  noting  that  the  expression  for 
the  longitudinal  velocity  v^,  Eq.  (IV-101),  is  of  the  Bame  fora  as 
Eq(lV-20),  the  potential  function  02n  and  ths  corresponding  longitudinal 
velocity  v^g  are  given  by  Eq.  (IV-22)  and  (IV-24)  respectively*  Using 
the  boundary  condition,  Eq.  (IV-25),  the  egressions  for  CQi  and  the 
velocities  vfc2  and  v02  are  given  by  Eq.  (!V-30)-(lV-32). 

To  satlofy  the  condition  of  zero  pressure  on  the  free  surface. 


the  fluid  potential  0^n  is-introduced.  The  expressions  for  0^n  and 


the  corresponding  velocity  components,  y^,  and  v^  are  given  by 
Eq.  (!V-33)-(lV-36)  where  the  coefficients  are  still  to  be  determined* 


The  coordinates  A^(t)  can  be  expressed  in  terms  of  q^t),  the  generalised 


coordinates  of  the  shell  motion.  This  is  done  by  utilisation  of  the 

Tr 

condition  that  the  fluid  surface  at  i  ■  ^  must  be  a  free  surface  with 


sero  pressure. 


To  determine  this  relation,  the  kinetic  energy  of  the  fluid, 


Eq.  (IV-3B)  is  evaluated: 


-103- 


Using  Eq.  (IV-106),  the  kinetic  energy  of  the  fluid  may  be  written  in 


terms  of  the  coordinate  q  (t)  : 


-forth  [  X  _ i _  +  - L - 

_  +*nhfaif&) 

'*  I  fur-mi] 


(jf-'ol) 


Substituting  Eq.  (IV-95)  for  into  Eq.  (IV-107),  the  kinetic  energy  of 
the  fluid  in  the  mode  "n"  is: 


Tn  ;/w  i  Tn  ,M  jfo  *  wjhfo 


where  the  coefficients  qn  and  t|n  are  respectively 


(w-ioa) 


n  =  -J-y  (3+^*)-7  ' 


\  #9  I _  —  j  _  -f&rt/i  fezzr)_h 

r  tn£  zrX**£  ~*1)  i^T 


•/-  2  y-/ 
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( jjr-'G?* ■) 


(&-m 


The  motion  of  the  fluid  surface  at  z  =  +  in  each  mode  n,  is 

required  for  the  analysis  of  the  response  of  the  partially  filled  tank 

to  dynamic  loading.  These  motions  can  be  obtained  by  an  integration 

with  respect  to  time  of  the  longitudinal  velocity  v  _  on  the  fluid 

z  j 

surface  z  =  +  Proceeding  as  in  Section  (IV-a),  the  fluid  die- 
placement  in  the  n  mode  (positive  upward),  ^  Q(r,  0,  t)  is  given  by 
the  relation 


CO 


As  in  the  case  of  the  three  constant  approximation,  ^  ^(r,  0,  t)  can  be 
written  as 

^ni<r»  6>  '  ’IT  **ni  Jn  <“lM  C0S  n0  (IV-112). 


Table  IV" 1 

Root1  for  J*'^) 


n  «  0 

B  *  1 

n  -  2 

1  -  1 

3.832 

1,841 

3.053 

1-2 

7.016 

5.332 

6.707 

1-3 

10.1? 

6.536 

9.97 

1  -  4 

13.32 

11.71 

13.17 

i  -  5 

16.47 

14.86 

16.31 

n  «  3 

n  -  4 

n  -  5 

n  ■  6 

4.20 

5.31 

6.40 

7.50 

7.89 

9.04 

10.52 

11.74 

11.17 

12.33 

13.99 

15.27 

14.37 

15.53 

17.27 

18.60 

17.52 

18.79 

20.53 

21.88 
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Cylindrical  Tanks* 

Expressions  for  obtaining  the  modes  and  the  frequencies  of  free 
vibrations  of  partially  full  cylindrical  tanks  are  presented  in  this 
section  of  the  report.  As  in  the  case  of  an  empty  tank,  two  sets  of 
approximate  displacements  are  considered]  l)  a  set  in  which  five  con¬ 
stants  are  retained  in  Eq.  (i-l)  to  (1-3)  and,  2)  a  set  in  which  three 
constants  are  retained.  While  the  latter  will  give  satisfactory  estimates 
of  frequencies  for  most  applications,  the  former  are  required  to  give 
accurate  values  of  strains  and  stresses  for  use  in  forced  vibration 
problems. 

Employing  the  Rayleigh-Hit*  method  as  in  Section  I,  the  total  kinetic 

▲u 

energy  of  the  partially  filled  shell  in  the  n  mode  is  given  by  the 
relation 


Tn  “  Tnshell  +  Tnfluid  'V_1' 

where  ^nBhell  is  given  by  Eq.  (1-7)  and  Tnfluid  is  given  by  Eq.  (IV-55) 
for  the  three  constant  approximation  or  Eq.  (IV-108)  for  the  five  con¬ 
stant  approximation. 

The  potential  energy  of  the  combined  system  is  not  effected  by  the 
fluid  in  the  shell  (the  effect  of  gravity  forces  on  the  potential  energy 
of  the  system  has  been  neglected  in  Section  (IV)  and  is  given  by  Eq.  (1-30)- 
(1-12). 

As  in  the  case  of  the  enpty  Bhell,  frequency  determinants  will  be 
given  for  the  case  in  which  both  membrane  and  bending  effects  are  considered 


and  for  the  case  of  a  membrane  shell.  For  the  mode  n  =  0,  special  condi¬ 
tions  prevail,  and  this  case  will  be  discussed  in  part  (f)  of  this  section. 


a)  Approximation  Using  Five  Constants  -  Membrane  and  Bending  Effects 


n  4  0 


Let  the  displacements  of  the  shell  be  given  by  Eq.  (l-4)-(l-6)  which 
are  repeated  here  for  convenience. 


u(z,  0,  t)  a  u(z,  0)elaJt  «  |u|  +  X(—  -  )  cosn  9 

l— a  a2  Ler  J 


(V*2) 


v(  2/  9 ,  t)  **  v(z,  *  V—  sinnd 

£1 


(v-3) 


w(z,  9,  t)  «  w(z,  a  Y  +  W(~  -  ~  )  coBn0 


(v-4) 


The  kinetic  energy  of  the  shell  mass  in  the  n  n  mode,  Tnslielj_>  is 


given  by 


2*  L 


Tnshell  “  2 


r  f  f  <“n  +  ^ +  »n> 0,1942 


(V-5) 


o  o 


and  the  kinetic  energy  of  the  fluid,  Tnfluid  is  given  by  Equation  (IV-108) 
for  the  five  constant  approximation: ^ 


nfluid 


m.  nacu  , 
~~7. - 


V3£2 

L  a. 


w2 +  v3  !r  ^ +  w3  S2  w  (v-6) 


(l)  For  free  vibrations,  the  time  dependency  of  the  displacements  is  taken 


icut 

as  e  . 
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Subst-1  <,uting  Eq.  (V-2)-(V~4)  into  Eq.  (V-5)  and  adding  Eq.  (V-6),  the 
total  kinetic  energy  of  the  partially  filled  shell  in  the  n  mode  is 
obtained  from  Eq.  ( V— l) : 

+  p  (1  +  ~T~  ^  >  **  *V*  ^  ™j  (T-7) 

vhere  the  coefficients  qn  ,  and  q^  are  given  by  Eqs.  (lV-109a)- 
(IV-109c)  respectively. 

The  potential  energy  Vn  stored  in  the  shell  can  be  expressed  as  a 
function  of  the  displacements  u(z,  6),  v(z,  0)  and  w(z,  0)  by  Eq.  (1-10)- 
(1-12) ,  which  upon  substitution  of  Eq.  (V-2)-(V-4)  become  the  expressions 
given  by  Eq.  ' ' -13)  and  Eq.  (I-l4). 

Applying  the  Rayleigh-Ritz  method,  the  following  set  of  five  homo¬ 
geneous  linear  equations  in  the  five  unknowns  U,  V,  W,  X  and  Y  are  obtained: 


nil  +  13V 


1  -V  n 


+  [f  [M  (1+  3  *  1  *  k  Y-'° 

(V-12) 

where  k  and  M  are  defined  by  Eq.  (1-21)  and  (l-22)  respectively,  7  is  the 

fractional  height  of  the  filling  &  ■  -  and  n  ,  ti  ,  and  n  are  defined 

a  'n  'n  n 

by  Eq.  (lV-109a)-(lV-109c) . 

Nonvanishing  solutions  of  Eq.  (V-8)-(V-12)  and  free  vibrations  exist 

only  if  the  determinant  of  the  system  vanishes.  This  leads  to  the 

determinantal  frequency  equation  shown  on  Page  96, 

For  given  values  of  n  and  7,  Eq.  (V-13)  yields  five  positive  roots 

M,,  defining  five  mutually  orthogonal  modes,  of  frequencies 
J 

_  M.Eh 

2  .  — J - -  J  -  1,  2,  3,  4,  5.  (V-14) 

J  ma2(l-i/2) 


J  -  1,  2,  3,  4,  5. 


(V-14) 


yL  m  height  of  filling 
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The  6hcpe  of  the  mode  pertaining  to  a  particular  frequency  u>.  can  be  found 

J 

U  V  X  Y 

by  computing  the  ratios  -  ,  ^  -  and  ^  from  any  four  of  the  Equations 

(V-8)  to  (V-12).  As  in  the  case  of  the  entity  tank,  in  general,  only  the 

lowest  frequency  cu^  will  be  required  and  for  simplicity,  the  subscript  J 

will  be  dropped  for  each  of  the  constants  U.,  V.,  W.,  X.  and  Y.. 

J  J  J  J  J 

Once  the  frequency  a>  and  the  ratios  of  the  constants  have  been  ev¬ 
aluated  for  a  particular  value  of  n,  the  displacements  of  the  shell  become: 


’  Cr.  [S  !  *  I  4  ■  0) 


cos  n8 


(Y-15) 


v(z,  6) 


w(z,  Q)  * 


C  -  - 

n  W  a 


0  2  + 

n  W 


sin  nO 


(i  _  j±  71 

'a  2a  'J 


cos  n 9 


(V-16) 


(V-17) 


The  constant  Cq  is  again  chosen  so  as  to  normalize  the  mode  to  the 
total  mass  of  the  empty  shell.  The  normalization  condition  thus  becomes: 


L  2  n 


mi  /  /  (un  +  vn  +  vn  >  ad0dz  +  2  *fluid  *  2naDni  “  Mi  (V’l8) 


where  T^uid  is  BPace  dependent  part  of  Eq.  (IV-108).  Substituting 
Eq.  (V-15)-(V-17)  into  Eq.  (V-18),  the  egression  for  the  normalization 


coefficient  C  is 
n 


C  » 
n 


•» 

L  1  -1/2 

(  J  (1,  .V.-8£!  )  (I)2  +  5  r3  £S  (1)1  l 

*  f  K  x  +  '  V  V  mi  V  j 


(v-19) 
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Using  Eq.  (V-19),  Eq.  (V-15)-(V-17)  give  the  normalized  modes  of  free 
vibrations  of  the  partially  filled  tank. 

The  normalized  strains  and  stresses  in  the  partially  filled  tank  ) 
may  be  evaluated  from  Eq.  (l-30)-(l-35);  using  the  ratios  of  the  constants 
and  the  expression  for  the  normalization  coefficient  that  have  been 
derived  in  the  present  section. 


b.  /'opyoxt^nti on  Using  Five  Consts.ntn-M<‘inbrane  Effects  Only. 

A.g  in  the  case  of  empty  shells,  a  membrane  approximation  may  be 
made  _n  which  the  strain  energy  of  bending,  V,,,  is  set  equal  to  zero. 

Fur  lain  shells  in  which  the  thickness  to  radius  ratio,  h/a,  is  very 
sty! i,  excellent  approximations  of  frequencies,  strains  and  stresses  can 
uc  "btained  from  such  a  procedure  for  modes  with  a  low  circumferential 
wave  number  "n"'.  It  may  be  noted  that  the  range  of  h/a  and  n  for  which 
t.b  ;y;;,brane  assumption  is  permissible  can  be  estimated  by  methods 
similar  to  those  of  Section  III  of  this  Report. 

The  frequency  equation  and  the  mode  shapes  for  the  partially 
IuIj.  membrane  shell  are  obtained  by  setting  the  value  of  the  constant 
L  equal  to  zero  in  Ec.  (V-7)-(V-ll) •  The  system  of  homogeneous  equation: 
then  becomes: 


+  l_p  {M(1  +  23L*^  '  1)J  Y  "  ° 

(V-24) 

The  determinantal  frequency  equation  is  obtained  by  setting  the 
determinant  of  the  above  system  equal  to  zero,  or  by  setting  k  °  0  in 
Eq.  (V-13)«  For  this  determinant  see  Page  101.  For  given  values  of  fc, 
n  and  y,  Eq.  (V-25)  yields  the  five  positive  roots  M. .  The  shape  of  the 
mode  pertaining  to  a  particular  can  he  found  by  computing  the  ratios 

-  and  ^  from  any  four  of  the  Eq.  (V-20)-(V-24) . 

Again  only  the  frequency  corresponding  to  the  lowest  root  M^,  will 
be  of  interest. 

Once  the  frequency  o>  and  the  ratios  of  the  constants  have  been  evalu¬ 
ated  for  a  particular  value  of  n,  Eq.  (V-l4)-(V-l8)  may  be  used  to  compute 
the  normalized  mode  shapes  of  the  tank.  The  normalized  strains  and  stresses 
in  the  partially  filled  tank  may  then  be  evaluated  from  Eq.  (l-30)-(l-35)> 
using  the  normalization  constant  CQ  and  the  mode  shapes  of  the  present 


section. 


ay)(3v-/) 
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Toxi ration  Using  Three  Constants  -  Membrane  and  Pending  Effects. 


n  ^  0 

Let  the  displacements  of  the  shell  be  given  by  Eq.  (l-42)-(l-44) 
which  are  repeated  here  for  convenience: 


u(z, 

0, 

t) 

-  u(z, 

9)  elmt 

„  Z  0  i<Ot 

a  U  —  cosne  e 
a 

(V-26) 

v(z, 

t) 

"  v(z, 

9)  e10* 

„  z  .  icot 

*»  V  -  sinnS  e 
a 

(V-27) 

w(z, 

t) 

-  w(z, 

9)  e1"* 

..  1  a  iCOt 

«  W  -  cosno  e 
a 

(V-28) 

The  kinetic  energy  of  the  partially  filled  shell  in  the  n “  mode  is 
given  by  Eq.  (IV-56)  where  m  is  given  by  Eq.  (IV-57)*  Substituting 
Eq.  (V-26)-(V-28)  with  Eq.  (V-5)  the  kinetic  energy  becomes 


m.nco  a  t3 


n 


a 


-  2^1  r  u2  ♦  v2  ♦  Ai  ♦  Ky ) 


(v-29) 


where 


m 


K 


vn 


n  m. 


«  22 

n  m^ 


(V-29a 


The  potential  energy  can  be  expressed  as  a  function  of  the  dis¬ 
placements  by  substituting  Eq.  (V-26)-(V-28)into  Eq.  (l-10)-(l-12),  thus 
obtaining  the  expressions  given  by  Eq.  (1-46)  and  (1-47). 

Applying  the  Rayleigh-Ritz  method,  the  following  set  of  three  homo¬ 
geneous  linear  equations  in  the  unknowns,  U,  V,  and  W  is  obtained. 


2M  -  4  -  n2(W)(l4x) 
L  r 


3n(3v,-l) 

V  + 

kn2(l-: 

26 

2 

(V-30) 

(l+3k) 

V  + 

2n  +  ^ 
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vbere  k  and  K  are  defined  by  Eq.  (l-2l)  and  (1-22),  and  7  is  the 

fractional  heickt  of  the  filling. 

Jtonvanishing  solutions  of  Eq.  (V-30)-(V-32)  exist  only  if  the 

detcnoinant  of  the  system  vanishes,  thus  leading  to  the  determinantal 

frequency  equation  shown  on  Page  lok. 

For  given  values  of  |,  n  and  7,  Eq.  (V-33)  yields  the  three  positive 

roots  M.,  defining  three  mutually  orthogonal  modes,  of  frequencies  given 
0 

by  Eq.  (V-l4) .  The  shape  of  the  mode  pertaining  to  a  particular  frequency 

U  V 

can  be  obtained  by  computing  the  ratios  ^  and  g  from  any  two  of 
Eq.  (V-30)-(V-32), 

As  in  the  case  of  the  empty  tank,  the  lowest  frequency  only  will 
be  required  and  the  subscript  J  will  be  dropped  for  simplicity. 


The  displacements  of  the 

shell  may  be  written  in  the  form 

u(z,  6)  m  C  H  —  cos.  n$ 
n  w  a 

(V-34) 

v(z,  G)  m  Cn  J  |  sin  n 9 

(V-35) 

w(z,  6)  m  ~  COS  tQ 

(V-36) 

where  the  constant  Cq  is  again  chosen  so  as  to  normalize  the  mode  to  the 

total  mass  of  the  empty  shell.  The  normalization  condition  is  given  by 

(2) 

Eq.  (V-l8),  which  upon  substitution  of  Eq.  (V-3^)-(V-36)  gives  the  follow¬ 
ing  ejpression  for  the  normalization  coefficient  Cq: 


(2)  For  the  three  constant  approximation,  Tfluld  is  the  space  dependent  part  of 
Eq.  (IV-55). 


i 


1+  Knf  +  $  + 


Using  Eq.  (V-37)>  Eq.  (V-3^)-(V-36)  give  the  normalised  modes  of  free 
vibration  of  the  partially  filled  tank. 

The  results  obtained  from  this  approximation  may  be  used  where  an 
e stir. ate  of  the  frequency  of  vibration  of  any  particular  mode  is  required 
However,  the  present  approximation  does  not  give  sufficiently  accurate 
results  for  the  strains  snd  stresses  in  the  shell,  and  the  approxima¬ 
tion  using  five  constants  must  be  employed  when  these  quantities  are 
required. 


d)  Approximation  Using  Three  Constants  -  Membrane  Effects  Only, 
n  4  0 

As  in  Part  (b)  of  this  section  a  membrane  approximation  may  be  used 
to  obtain  the  frequencies  of  sufficiently  thin  shells  in  the  lover 
modes  of  n.  The  frequency  equation  and  the  mode  shapes  for  the  membrane 
shell  are  obtained  by  setting  the  value  of  the  constant  k  equal  to  tero 
in  Eq.  (V-30)-(V-33)«  The  determinantal  frequency  is  given  on  Page  107 . 

The  remarks  of  Part  (c)  of  this  section  regarding  the  validity 
and  applicability  of  the  three  constant  approximation  also  hold 


for  this  case. 


7L  -  height  of  filling 
(V-33a) 
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e’)  Determination  of  Approximate  Frequencies  of  Partially  F»n 
Cylindrical  Tanks  -  Rayleigh’s  Method . 


If  the  frequencies  and  mode  shapes  of  on  empty  cylindrical  tank  are 
known,  an  expression  for  the  determination  of  approximate  values  of  the 
frequencies  of  the  partially  full  tank,  in  terms  of  these  quantities 
can  he  derived  from  Rayleigh's  principle.  For  a  partially  full  tank, 
Rayleigh's  principle  may  he  written  for  the  n  n  mode  of  vibration: 


T 


npax 


rp  ^  *1* 

“nshell,  max  nfluid,  max 


V 

n,  max 


(V-38) 


where  the  subscript  "max"  indicates  the  maximum  value  of  the  respective 
quantity  and  the  Tncj)ell^  Tnfiuid  and  vn  are  given  by  Eq.  (1-7),  (IV-55) 
and  (1-10)  respectively.  Writing  T  in  the  form 


T  «  a)2  T 

n,  max  n  n,  max 


(V-39) 


where  Tq  is  obtained  by  substituting  the  displacements  u(z,  0),  v(*,  6) 
and  v(z,  6)  instead  of  the  velocities,  into  the  equations  for  kinetic 
energy.  substituting  Eq.  (V-39)  into  Eq.  (V-38),  the  frequency  of 
the  n  mode  of  vibration  of  the  partially  full  shell  can  be  written 
in  the  following  form: 


0) 


V 

n,max 

V 

nshell,  max 


$ 

_ nshell,  max _ 

nBhell,  max  +  ^nfluid,  max 


(V-40) 


Substituting  the  mode  shapes  for  the  empty  shell  into  the  expressions 

^0r  ^n  ^n  aad  noting  that  for  this  case,  the  frequency 

of  the  n^  mode  of  the  empty  tank,  a>ne,  is  given  by 
w 

2  n.max 

03  m  — — Z-tzr. 

no 


T 

nshell, max 


(V-4l) 


the  frequency  equation,  Eq.  (V-i+O)  becomes 


CD 


;  2 

O  CD 

n  ne 


(V-42) 


where 


nshell,  max 


n 


T  +  T 

nshell,  max  nfluid,  max 


(y-43) 


The  value  of  &n  for  the  five  constant  approximation  is  given  by  the 


relation. 


6  = 
n 


c2sV 


1  +  -A -  7  fiS  + «  62(A)2  +  <?  62  )U  (v-44) 

6  [  (n  fn  m^Hr  ( n  ra^  W  J  '  ' 

where  Cn  and  n,  rjf^  are  given  by  Eq.  (1-29)  and  Eq.  (lV-109a)-(l09c) 
respectively. 

The  value  of  6^  for  the  three  constant  approximation  is  given  by 


the  relation 
6 


n 


C2K  t273 
n  n* 


(V-45) 


1  + 


where  Cn  and  Kq  are  given  by  Eq.  (1-56)  and  Eq.  (V-29a)  respectively. 

Equations  (V-42)-(V-45)  may  be  used  to  obtain  approximate  values  of 
the  frequency  <an  of  the  partially  filled  tank  for  modes  in  which  n  >  2. 

It  should  be  emphasized  however  that  for  the  mode  n  ■  1  and  for  all  modes 
in  which  the  strains  and  stresses  in  the  shell  are  requir  ed,  the  value  of 
the  frequency  should  be  computed  from  the  determinantal  frequency  equations 
given  in  Sections  (V-a)-(V-d)  and  the  correct  mode  shapes  should  be  evalu¬ 
ated  from  the  corresponding  equations  given  in  these  sections. 


Table  (V-l)  shows  the  application  of  Equations  (V-4l)«(V-45)  to 

the  cases  of  9/10  full  steel  and  concrete  tanks*  The  value  of  the 

o 

frequency  number  M  ,  which  is  proportional  to  u>  is  evaluated  from 
n  n 

the  relation 


M  »  8  M 
n  n  ne 


(V-46) 


where  Mng  is  the  frequency  number  for  the  empty  shell.  It  is  seen 
that  good  approximatiorB  to  the  frequency  are  obtained  in  both  cases. 
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f )  n  »  0 


The  frequency  <dq  of  *  partially  full  tank  can  be  determined  in 
terms  of  the  frequency  to  of  the  empty  tank.  Noting  that  the  displace- 
ment  of  the  shell  in  the  mode  n  «  0  is  purely  radial,  the  total  effect 
of  the  fluid  filling  is  an  increase  in  the  mass  of  the  shell,  m^. 

Considering  the  s-lndependent  radial  motion  of  the  shell  as 
described  in  Section  (lV»b), 

v  -  qo(t)  (V-V?) 

and  using  Eq.  (IV-88)  and  (1 V-89),  the  Increase  in  the  mass  of  the  shell 
due  to  the  fluid  filling  of  height  fL  is  given  by  The  frequency 

of  the  partially  full  fuel  tank  can  then  be  computed  from  the  relation 


vhere 


(V-49) 


Tnble  V-l 


V  - 

-S- 


Determination  of  M  «•  &  M 
n  nte 


Unprotected  Steel  Tank  9/10  Full 


L/a  -  0.8,  |  « 


1 

1200 


7  *  *9 


n 

M 

ne 

Eq.(l-lfl) 

B 

n 

Eq.  (V-44) 

M 

n 

Eq. (V-46) 

Mn 

Eq.  (V-25) 

2 

.1906 

.0495 

.00943 

.00809 

3 

.0859 

.0578 

.00497 

.00475 

If 

.0432 

.0659 

.00285 

.00281 

5 

.0241 

.0741 

.00179 

.00175 

6 

.0147 

.0816 

.00120 

.00120 

Protected  Steel 

Tank  9/10  Full 

L/a  -  0.8,  |  - 

31  •  7  •  -9 

n 

Mne 

Eq. (1-23) 

6 

n 

Eq. (V-44) 

Mn 

Eq.(V-46) 

M 

n 

Eq.  (v-13) 

2 

.1911 

*3859 

•0737 

.0688 

3 

.093^ 

.4229 

.0395 

.0387 

If 

.0659 

.4560 

.0301 

.0301 

^  V  E*.  v ->>v  -r -7><**~*  r****: *• -ft. 
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VI.  Modes  and  Frequencies  of  Cylindrical  Tanks  with  Windgirders 

In  Section  (I)-Section  (V)  of  this  report,  approximate  expressions 

for  the  frequencies  and  modes  of  empty  and  partially  filled  cylindrical 

storage  tanks  were  developed.  In  obtaining  these  expressions,  the  effect 

of  the  wind  stiffening  ring,  i.  e.  windgirder,  which  is  generally  placed 

at  or  near  the  top  of  the  structure  was  neglected.  The  effect  of  the 

windgirder  on  the  frequencies  obtained  in  the  previous  Sections  can  be 

derived  from  a  procedure  using  Rayleigh's  Method.  An  upper  bound  to  the 

frequency  of  the  tank  with  a  windgirder  in  terms  of  the  frequency  o> 

of  the  storage  tank  without  a  windgirder  is  obtained.  It  will  be  shown 

■that  for  the  type  of  tanks  under  consideration,  about  40  feet  high-100  feet 

diameter,  the  effect  of  the  windgirder  on  the  tank  frequency  is  very 

small.  In  general,  it  may  be  neglected  in  computations  and  the  formulas 

presented  in  the  previous  Sections  of  this  report  can  thus  be  used  for 

the  determination  of  frequencies  and  modes  of  free  vibrations  of  struc- 

(2) 

tures  with  windgirders .  ' 


(l)  The  dimensions  and  the  type  of  windgirders  in  the  analysis  are  taken 
from  the  American  Petroleum  Institute  Specification  for  Welded  Oil  Storage 
Tanks,  American  Petroleum  Institute  Standard  12C,  Fifteenth  Edition, 

March  1958,  American  Petroleum  Institute,  New  York,  Pg.  13  and  70. 


(2)  ThiB  does  not  imply  that  the  blast  loading  on  structures  with  wind- 
girders  may  not  differ  from  the  loading  without  ouch  a  girder.  Information 
on  blast  loading  of  tanks  with  windgirders  is  not  yet  available. 


Consider  a  storage  tank  with  a  windgirder  attached  to  the  tank  at  a 


height  L  from  the  tank  bottom.  Fig.  VI-1.  If  the  frequencies  and  mode 

shapes  of  an  empty  or  partially  full  cylindrical  tank  without  a  windgirder 

are  known,  an  expression  for  the  frequencies  of  the  stiffened  tank  is 

derived  as  follows.  Rayleigh's  Principle,  T  »  V  .  can  be  written 

^  max  max' 

for  the  nta  mode  of  vibration  in  the  form: 


Tns  +  TnR 


V  +  V  n 
ns  nR 


(VI-1) 


where  T  and  V  are  the  maximum  values  of  the  kinetic  energy  and 
ns  ns 

potential  energy  respectively  of  the  empty  or  partially  filled  unstiffened 
shell  and  and  V^  are  the  maximum  values  of  the  kinetic  energy  and 
potential  energy  of  the  windgirder.  Defining  the  quantity  Tq  as  in 
Eq.  (V-39),  the  frequency  of  a  shell  with  a  windgirder  is  obtained  from 
the  relation: 


a 3 


ns 
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+  V 


nR 


ns 


+  T 


nR 


(VI-2) 


—  1  12 
By  definition,  T  „  »  x  M  „  and  therefore  V  ■  *  M  co  .  where  M  „  and 
'  ns  2  ns  ns  d  ns  ns  ns 

till 

o>n0  are  the  generalized  mass  in  the  n  mode  and  the  frequency  in  the  n 
mode  respectively  of  a  shell  without  a  windgirder.  Substituting  the  above 
expressions  into  Eq.  (VI-2),  the  frequency  is  given  by: 
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where  M  0  is  the  generalized  mass  of  the  windgirder. 
rut 


(VI-3) 


Assuming  that  the  shapes  of  the  modes  of  the  tank  with  the  wind- 
girder  are  the  same  as  the  shapes  of  the  unstiffened  tank,  the  windgirder 
ring  must  undergo  the  displacements  w(L,  0 )  and  v(L,  0),  of  the  tank  at 
the  height  z  «  Lj 


v(L,  0)  *>  C 


ris  n 

n  L.VJ  a  J 


sin  n0 


(VI-4) 


w(L,  0)  n  Cn  £  +  (^T1)  I  cos  n6  (Vl-5) 


$ 


where  C  ,  the  normalization  coefficient  for  the  unstiffened  tank  and  the 
n  V  Y 

mode  shape  ratios  ~  and  ~  are  given  in  Section  (I)  or  Section  (V). 

n  n 

The  potential  energy  of  the  ring,  V  in  terms  of  the  ring  displacements 
v  and  w  iaj 


Vr 
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f  (vee  *  w>i 
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d0  +  f  (v0  -  w)2  d0 


o  o 

which,  upon  substitution  of  Eq.  (VX-4)-(VI-5)  becomes 
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(VI-7) 


The  generalized  mass  of  the  ring,  M^,  is 


which,  upon  substitution  of  Eq.  (VI-4)-(VI-5)  becomes: 

_2^*,VnLx2  fYn  ,?L  -  L.  I  ^ 


a  „  «2  ( /n  Lx2  [~Yn 

,pRVaCn  +  r 

(  n  —  n 


/PL  -  LN 
+  (  575  ) 
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Since  in  all  cases,  the  modes  of  vibration  of  the  un stiffened  shell  have 
been  normalized  to  =  2ital/n^,  the  total  mass  of  the  empty  tank,  the 
generalized  mass  of  the  shell  in  all  modes  "n"  is  given  by: 


M  „  =  2rtaLm. 
ns  i 


(VI-10) 


Once  the  frequency  of  the  unstiffened  shell,  co  ,  is  known,  Eq.  (VI-7)- 

ns 

2 

(VI-10)  can  be  substituted  into  Eq.  (VI-3)  and  the  value  of  to  can  be 
computed. 

Table  (VI-l)  contains  the  values  of  the  frequencies  u>n  and  o>Qa  for 

the  typical  case  of  an  empty  steel  shell,  100  ft.  in  diameter  and  40  ft 

high,  having  the  windgirder  section  shown  in  Fig.  (VI-2).  The  vindgirder 

has  been  attached  to  the  shell  at  a  height  L  -  39  ft.  It  is  seen  that  the 

correction  due  to  the  windgirder  is  very  small;  this  correction  is  less 

than  5  1/2  per  cent  in  all  modes.  The  situation  in  fluid  filled  tanks  is 

quite  Bimilar.  Therefore,  it  is  concluded  that  for  the  type  of  storage 

tanks  under  consideration,  the  effect  of  the  windgirder  may  be  neglected 

and  the  formulas  of  Section  (i)-(v)  can  be  used  even  for  structures  with 
(4) 

windgirders. 


[3) American  Petroleum  Institute  Specification  for  Welded  Oil  Storage  Tanks 


Fifteenth  Edition,  March  1958/  Pg.  13  and  70:  The  minimum  required  section 

p 

modulus  for  a  stiffening  ring  is  Z«  0.0001  D  Eg  where  2  is  the  section 
modulus  in  in  ,  D  is  the  nominal  diameter  of  the  tank  in  feet,  and  is 
the  total  height  of  the  tank  shell  in  feet.  For  a  tank  of  100  feet  diameter 
and  40  feet  height,  Z  a  40  in3.  The  section  modulus  of  the  windgirder  shown 
in  Fig.  (VI-2)  is  47.7  in3  (See  Pg.  70,  A.P.I.  Specs). 


(4)  See  Footnote  2,  Section  VI 


Table  VI-1 


Effect  of  Wlndglrder  on  the  Frequency  of  Cylindrical  Tanka 

Empty  Steel  Tank 

—  p  |i 

L  -  40  ft«,  a  -  50  ft.,  L  «  39  ft.,  »  11.2  in  ,  IR  a  449  in 


a 

“ns  (rad/sec) 

o>n  (rad./sec) 

#  Difference 

1 

217.4 

211.5 

2.7# 

2 

140.3 

134.2 

4.4# 

3 

94.2 

89.4 

5.1# 

4 

66.8 

63.4 

5.2# 

5 

49.9 

47.9 

4.1# 

V  *• 

■A.  ~L 
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» * '  ■ '  v- » • . 
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Forced  Vibrations  of  F  ntv  end  Partially  Full  Cylindrical  Fuel  Tanks. 


The  respo.iccsof  empty  and  partially  (up  to  9/10  of  the  height)  full 
cylindrical  tanks  to  dynamic  loading  are  analyzed  in  this  Section.  For 
empty  shells,  the  blast  loading  consists  of  a  radial  pressure  dynamically 
applied  to  the  shell,  while  for  partially  full  shells,  the  vertical  pres¬ 
sures  of  the  blast  on  the  surface  of  the  liquid  must  also  be  considered. 
Information  regarding  the  pressure  history  is  given  in  the  various  Armour 
Institute^Pcports  and  specific  methods  for  the  utilization  of  the  pressure 
data  are  developed  in  the  present  paper.  The  cases  of  a  partially  full 
tonic  and  an  empty  tank  will  be  separately  treated. 

Using  the  modes  of  free  vibrations  of  the  empty  or  partially  full 
tank  respectively  as  generalized  coordinates,  the  response  can  be  expressed 
in  terms  of  the  infinite  number  of  modes.  The  required  frequencies  and 
modes  of  empty  and  partially  full  cylindrical  shells  were  determined  in 
Sections  (I)  and  (V)  of  the  report  respectively.  For  each  integral  value 
"n"  of  circumferential  waves  in  the  shell  displacement,  there  exists  an 
infinite  number  of  frequencies  J  -  1,  2,  3*  •  •  •  and  corresponding 
mode  shapes  u^z,  0),  v^(z,  0)  and  ^(z,  0)»  Calling  qQj(t)  the 
generalized  coordinate  corresponding  to  the  mode  "nj",  the  response  of  the 
shell  can  be  written  in  terms  of  a  summation  of  the  normal  modes: 

00  CO 

u(z,  0,  t)  -  £  £  UqJ(z,  SJq^t)  (VII  -1) 

DaO  J“1 


(l)  Blast  Effects  on  Storage  Tank  Type  Structures,  Armour  Research 
Foundation  Report  Number 


(vir  -2) 


v(z,  e,  t)  =  2^  ^ 
n=0  j=l 

00  CO 

w(z,  e,  t)  -  £  £  vnj(z,  ejq^Ct)  (VII  -3) 

ifco  j=i 


Yflien  determining  the  frequency  numbers  (which  are  proportional 

O 

to  a>nj  )  it  is  found  that  the  values  are  smaller-usually  very  much 
smaller,  than  unity,  while  the  values  of  for  J  >  1  are  larger  than 
unity.  The  modes  characterized  by  j  >  1  represent  primarily  extensional 
high  frequency  oscillations  of  the  shell.  For  loads  of  moderately  long 
duration  compared  to  the  fundamental  period  T  .  for  each  value  of  n, 
these  high  frequency  modes  give  only  small  contributions  to  the  response 
of  the  shell.  The  major  response  of  the  structure  can  be  determined  from 
the  modes  characterized  by  J  «  1  and  all  modes  for  which  J  >  1  are  there¬ 
fore  dropped  from  the  summation. 

The  special  case  n  a  0  has  been  discussed  in  Sections  (i-e)  and(V-f) 

1 

of  the  report.  For  the  empty  shell,  an  infinite  number  of  frequency  numbers 
MoJ  are  found  to  lie  between  closely  spaced  limits,  thus  making  a  model 
analysis  for  the  case  n  =  0  impracticable.  A  suitable  alternative  method 
to  obtain  the  contribution  of  the  n  =  0  term  is  to  consider  the  shell  acting 
as  a  series  of  separate  rings,  each  having  a  purely  radial  displacement 
vQ(t).  For  the  fluid  filled  shell,  however,  the  displacements  of  the  rings 
on  different  levelB  are  coupled  by  the  fluid  and  a  spectrum  of  gradually 
increasing  frequencies  is  again  obtained.  The  fundamental  frequency  J  n  1 


in  the  filled  shell  contains  again  the  major  contribution  to  the  response 
for  n  o  0,  as  in  the  case  where  n  /  0. 

The  analysis  of  the  response  of  the  shell  in  the  mode  n  «  0  will 
be  discussed  separately  in  Part  (  C  )  of  the  present  section. 

For  the  pressure  loadings  of  moderately  long  duration  under  con¬ 
sideration,  the  response  of  the  shell  can  therefore  be  written  in  the  form 

00 

u(z,  8,  t)  a  V  UQ(z,  0)qn(t)  (VII  -4) 

n=rl 

00 

v(z,  e,  t)  =  V  vn(z,0)qn(t)  (vir-5) 

n=l 

00 

w(z,  6,  t)  =  Wo(z)qQ(t)  +  y  Vn(z,  ejq^t)  (VII-6) 

n=l 

where  the  functions  u  ,  V  and  V  are  the  modes  of  free  vibrations  of 

n  n  n 

the  empty  or  the  partially  full  shell  corresponding  to  J  »  1.  Both  the 
five  constant  and  the  three  constant  approximations  of  Section  I  will  be 
considered  in  the  following  analysis.  It  should  be  noted  that  for  cases 
in  which  shell  strains  and  stresses  are  required,  the  five  constant  analy¬ 
sis  should  be  used.  Expressions  will  now  be  developed  for  the  determina¬ 
tion  of  the  generalized  coordinates  ^(t)  from  the  information  given  in 
the  Armour  Report. 


a)  Plast  Pressures  on  Fuel  Tanks  -  Nomenclature  and  Measuring  Arrangement 


The  dynamic  load  information  required  to  determine  the  response  of 
a  cylindrical  fuel  tank  under  blast  loading  is  obtained  from  a  report  on 
shock  tube  experiments  by  the  Armour  Research  Foundation  (See  footnote  1, 
Section  7)*  A  description  of  the  shock  tube  tests  is  given  in  the  Armour 
report.  The  gages  were  arranged  to  give  enough  pressure  information  to 
allow  for  a  meaningful  forced  vibration  analysis  of  the  tank.  Two 
types  of  pressures  were  required;  a)  the  radial  pressures  on  the  tank 
walls,  and  b)  the  vertical  pressure  on  the  surface  of  the  fluid  filling. 
In  this  section,  the  nomenclature  and  the  arrangement  of  the  pressure 
information  required  for  a  forced  vibration  analysis  will  be  described. 

For  purpose  of  obtaining  vertical  pressures  on  the  fluid  surface  of 
the  tank,  the  surface  is  divided  into  sixteen  sectors  by  means  of  eight 
diametral  planes  spaced  22  l/2  degrees  apart.  Each  sector  is  denoted 
as  S  ,  k  =  1-16  and  is  located  by  an  angle  0  measured  to  the  midpoint 

K 

of  sector  k.  The  intersection  of  the  bisector  of  8^  and  the  periphery  of 
the  shell  is  called  i,  1  0-15.  The  angle  0  is  measured  positive 

clockwise  from  the  zero  position  which  is  the  point  at  which  the  shock 


wave  first  hits  the  shell.  The  sector  8^  is  bounded  by  the  angles  0^  and 


Each  sector  is  further  divided  radially  into  three  subareas  8^,  where 


the  subscripts  p  **  1,  2,  3  refer  to  the  outermost,  middle  and  interior 

subareas  of  the  sector,  respectively.  Each  subarea  contains  pressure 

gages  from  which  a  mean  vertical  pressure  is  obtained.  The  pressure 

P,  is  considered  constant  over  the  subarea  8,  .  Fig.  (VII-2) 
kp  kp 


and  (VII-3)  show  the  divisions  of  the  fluid  surface  and  the  nomenclature 
used  in  the  following  sections  of  the  report. 

For  purposes  of  obtaining  radial  pressures  on  the  walls  of  the  tank, 
gages  were  placed  to  measure  the  average  pressure  along  a  generator 
denoted  by  i  at  the  various  heights  shown  in  Fig.  (VII-l) .  The  outside 
wall  pressure  are  thus  obtained  for  a  sector  in  terms  of  the  four 
pressures  p  ^-p  ^or  pressure  on  the  inner  walls  of  the 

shell,  the  pressure  P^  ia  used  for  a  9/10  full  tank  and  the  pressures 
P^  and  P^  are  used  for  a  half  full  tank.  In  each  subarea  shown  in 
Fig.  (VII-l),  the  pressure  p  ,  h  =  1-6,  is  consiuered  constant  over  the 
area. 

Once  the  weapon  yield  and  the  peak  ambient  overpressure  are  selected, 
the  external  pressures  acting  on  the  tank  can  be  obtained,  and  the 
forced  vibration  analysis  of  the  following  sections  can  be  made. 
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b)  n  ^  0  Five  Constant  Approximation  ~  Partial  Filling  0  <  y  <  1.  (See  Flg»  VII-f 


Substituting  the  expression  Eq.  (V-15)  to  (V-17)  for  the  mode  shapes 
into  Eq.  (VII-4)-(VII-6)  for  the  response,  the  displacement  u,  v,  v  of  a 
point  on  the  shell  becomes: 


u(*,  e,  t) .  f  c  r  l +  r  (4  •  )  °°m$  ,B(t)  <vn-7) 

n  n  a  4a 


"  y 

v(z,  6,  t)  »  \  C  ~  |  sinne  q  (t) 
i—i.  n 


(VII-8) 


»  Ty 

v(z,  0,  t)  .  Vo(»)q0(t)  +  V  CQ  +  (£  -  ~  )  cosne  qn(t).  (VII-9) 

L  n  J 


Each  mode  is  normalized  to  the  total  mass  of  the  empty  tank,  2xaLr&i 
and  the  normalization  coefficients  Cq,  are  given  by  Eq.  (V-19). 


The  generalized  coordinates  qQ(t)  satisfy  Lagrange's  equations 


Si*  °n  ^n  *  5 


2naLm. 


(VII-10) 


where  o)q  ,  the  frequency  of  vibrations  of  the  tank  in  the  tl  mode  is 
computed  from  the  lowest  root  of  the  frequency  determinant  Eq.  (V-13)> 

Qor  Eq.  (V- 25)  if  bending  effects  are  small  such  that  the  shell  can  be 

considered  to  act  as  a  membrane]]and  Q^t)  is  the  generalized  force  in  the 

th  , 
n  mode. 

The  generalized  force,  Q^t),  may  be  written  as  the  sum  of  two  terms: 

^(t)  -  Q^t)  ♦  tyt) 

where  Q^t)  is  the  generalized  force  produced  by  the  radial  blast 


(VII-11) 
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pressure  on  the  tank  walls  and  4^_(t)  is  the  generalized  force  produced 
by  the  vertical  blast  pressures  acting  on  the  surface  of  the  tank  filling. 
The  force  Q^(t)  is  given  by 
2*  L 

Q^(t)  =  J'  J  P(0,  z,  t)  vn(0,  z)  ad0dz  .  (VII-12) 


0  o 


Substituting 


w  (z,  6)  =>  C 
n'  ’  7  n 


z  /  n 

a  +  W 
_  n 


—  ) 
2a  ' 


cosn0 


(VII-13) 


into  Eq.  (VII-12),  Q^(t)  may  be  written  in  the  form 


2* 

■/  , 

r  L 

f  He,  *,  t)  cn 

i  +  (!s - h\ 

a  'W  2a' 

dz  f 

i  1 

J 

L  0 

L  n  _J 

J 

Referring  to  Section  (Vll-a)  in  which  the  pressure  measurements  made  on 
the  Armour  Shock  Tube  model  are  described  in  detail,  it  is  evident  that 
the  spatial  integration  with  respect  to  z  can  be  performed  at  each  of  the 
sixteen  vertical  generators  "i"  on  the  shell,  i  -  0,  1,  2,  .  .  .  16  as 
shown  in  Figure  (VII-l). 

Denoting  the  integral  for  each  i  as  f^, 

Y 


f„i  "  /  Pl<*’  Cn  !  + 

w  —  n  _ 


dz 


(vil-15) 


and  letting  P  ^  be  the  pressure  at  the  various  heights  h  of  the  tank 
(h  a  1,  2,  3 >  5 ,  6)  in  the  segment  8  containing  the  point  i,  the  func¬ 

tion  ffii  is  easily  computed  in  terms  of  the  P^.  For  example,  for  the 

(2) 

case  of  a  9/10  full  tank,  7  ■  0.9,  f^'  becomes 


f  ,  -  C 
ni 


na{2  Q  .105PU  -  .015P12  +  .075P13  +  .045(PlU  -  Pl67J 

_pll  +  P12  +  P13  *  °-333  (p4fc  -  Pt6)] 


♦V  -31  5* 


(YII-14) 


(VII-16) 


Equation  (VII-lU)  may  be  written  in  the  form: 

2n 


^(t)  =  JQ  fQ(0,  t)  cosn0  ad 0 


(VII- 


where  fn(0,  t)  is  a  continuous  function  of  the  angle  0.  The  function 


fQ(0,  t)  ia  expressed  as  a  Fourier  Cosine  Series, 


8 


f  ( 9 ,  t)  «*  a  (t)  + 
n'  no 


1=1 


anj(t)  COBJ0 


(VII- 


where  the  coefficients  a  and  a  .  are  determined  so  that  f  (0.  t)  *  f  . 

no  nj  n'  7  '  ni 

when  0  «=  0£  (i  a  1,  2,  ...1 6.)  This  requirement  restricts  the  series  to 


nine  terms  and  the  coefficients  a 


nj 


(3) 


are 


a  “IS 


no 


(VII. 


i=>0 


15 

an,J  c  B  £fni(if)cos^ 
i=0 

j  ■  1,  2,  3/ ' 


(VII- 


an8  "  IB  fni(i!)  COBnL 


(VII. 


i=0 


( 2  )  For  the  case  of  a  5/10  full  tank,  the  formula  for  f  ^  becomes 

fnl  ”  Cnal2  [_-  -105PU  '  -015I,i2  +  •OT5p13  +-«.i*5(PlV  Pl6)-  -0BP, 


+  Cn.a  W 


n 


+  P12  ^  ?i3  t  0.333(2^  “  ^ig)  "  1‘333P^ 


(  3 )  The  harmonic  analysis  in  which  a  Fourier  Series  is  passed  through 

the  sixteen  points  i  and  where  the  coefficients  a  .  are  determined  is 

nj 

given  in  the  Appendix  to  Bection  VII. 
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Sub  otituting  Eq.  (VIII-lQ)  into  Eq.  (VII-17), 

8 

+ 


o+  I  a^tjcosje  coan0  ad0  (VTI-22) 


no 

o  J=1 


and  uaing  the  orthogenality  relation 
2  n 


/ 


cosj0  cosn0  d0 


0  n  j  j 


n  n  »  j 


(yii-23) 


the  generalized  force  Q^t)  becomes 


V1’ 


«  Jta  a  (t)  « 

nn'  ' 


n  *  1 y  2,  • • .8 


(VII-24) 


The  load  information  given  in  the  Armour  report  is  sufficient  to  give 
generalized  forces  in  the  nine  modes,  n  ■  0  to  n  *  8.  In  general  there¬ 
fore,  a  nine  mode  elastic  analysis  may  be  made  for  the  determination  of 
the  response  of  the  shellj  in  practice,  vith  loads  of  moderately  long 
duration  being  considered,  the  series  can  usually  be  stopped  after  the 
fourth  to  sixth  mode. 


The  generalized  force  Q^t)  due  to  the  vertical  pressure  acting  on 

(4) 


the  fluid  surface  at  z  »  yL  is  given  by 

2  »t  a 


(^(t)  »  -J  J  P(r,  0,  t)jn(r,0)rdrd0 


(VII-25) 


where J-n(r,e)  *s  the  space  dependent  part  of^.n(r,  6,  t)  which  is  the 
displacement  of  the  fluid  surface  at  z  o  yh.  From  Eq.  (IV-110)-(IV-112) t 


*.<*•>- X  l 

V  1-1  °  i«l 


tt.if 


COSD0  (VII-26) 


(4)  The  minus  sign  is  due  to  the  fact  that  the  pressure  P(r,0,  t)  is 
considered  positive  when  acting  downward,  while  positive  displacements 
J^n(r,  0)  are  up  weirds. 


j, 
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vhere 


a  .  ?L 

l  2anitanh  (~-)  1  , 

'  nl  n'  ni'  L_  a 


aH7L 
tanh(-~ -  ) 


+  tanh( 


V’’1,  \ 

2a  •) 


a  .7t 
8lnh(-~ — ) 


( VI I -27 ) 


(Substituting  Eq.  (VII-26)  into  Eq.  (VII-25),  {^(t)  becomes 

2it  a  » 

Q^Ct)  55  “  /  /  p(r>  0>  t>  Cn  2t[Z^niJn(”F~)1  0080  rdrd£?  <m“28> 

o  O  -> 

Referring  again  to  Fig.  (VII-2)  and  (VII-3)  of  Section  (Vll-a),  the 
fluid  surface  is  divided  into  circular  sectors  denoted  by  k,  k  *  1,  2,  3; • • «l6* 
Moreover,  each  sector  k  ia  further  subdivided  into  the  subareas  3^,  p  =  1, 

2,  3,  corresponding  to  the  location  of  the  roof  preBBure  measuring  gages. 


In  each  subarea 


k  ■  1,  2  f  ...  16 

p  =  1,  2,  3 


(VII-29) 


the  pressure  is  considered  constant  and  the  integration  of  Eq.  (VII-28) 
can  be  performed  analytically  in  each  subarea.  The  expression  for  the 


generalized  force  ^(t)  becomes  therefore 
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The  contribution  to  Q^t)  from  the  vertical  pressures  acting  on  the 


subarea  in  Eq.  (VII-30)  can  be  written  in  the  form 
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where  the  limits  on  the  integrals  are  determined  by  the  particular  subarea 


S  .  Interchanging  the  integration  with  respect  to  r  and  the  summation 


over  i,  and  defining  the  quantities 
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(VII-33) 
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Analytical  expressions  for  the  integrals  1^ 
Table  (VII-l)  for  n  -  0(1)6. 

The  generalized  force  Q^t)  is  given  by 
16  3 

E%(t) 


[]Eq.  (VH-3^ 


are  given  in 


(VII-35) 
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which,  upon  substitution  of  Eq.  (VH-3^)  becomes: 
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In  practice,  a  sufficient  number  of  the  fluid  modes,  denoted  by  the  sub¬ 
script  "i"  must  be  considered,  commensurate  with  the  required  accuracy. 


For  the  type  of  study  considered,  the  first  three  terms  of  the  series  of 
fluid  modes  suffice  and  the  upper  limit  in  the  summation  over  i  in 
Eq.  (VII-36)  becomes  three. 

The  generalized  force  ^(t)  can  be  evaluated  from  Eq.  (VII-ll)  and 
the  generalized  coordinates  q^t)  can  be  computed  by  a  numerical  integra¬ 
tion  of  Eq.  (VII-10).  A  suitable  numerical  integration  technique,  due  to 
Noumcrov, ^  evaluates  q^  by  a  forward  step  integration  in  time.  Writing 
Eq.  ( VII-10)  and  the  initial  conditions  for  a  system  starting  from  rest, 
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(VII 


qn(0)  o  q^O)  =  0  ' 


Let  "k"  be  the  interval  of  the  time  steps,  the  recurrence  formula  for 
qfl(t+k)  becomes 
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(VII 


Equation  (VII-38)  allows  the  determination  of  the  generalized  coordin 
ate  q^t+k)  in  terms  of  known  values  at  two  previous  time  steps.  The 
integration  is  started,  using  the  formulas 
qn(0)  »  0 


(5)  Numerical  Methods  in  Engineering  by  M.  0.  Salvadori  and  M.  L.  Baron, 
Prentice  Hall,  Second  Printing,  1955*  Pg»  H8  ff 
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-  ^F(0) 


2  2 
+  2k  a) 
n 


(vn-39) 


Once  the  generalized  coordinates  q^t)  have  been  evaluated,  the 


response  of  the  shell  is  obtained  from  Eq.  (VII-  ^  )  (VII-  9  ).  Using 

Eq.  (l-30)-(l-32),  the  strains  in  the  shell  are  given  by  the  expressions^ 
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Using  Eq.  (VII-4o)-( VII-42),  the  membrane  shell  stresses  can  be  com¬ 
puted  from  the  relations: 
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°z  0  =  aH+y)  _z0_ 


(VII-45) 


(6)  The  strain  eQg  in  the  mode  n  «  0  is  given  by  Eq.  (VII-  73  )  in  Part  (c) 


of  this  section. 


V-V-V  /• . 
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For  the  c.'t.'io  of  stool  shells  with  concrete  shielding,  bending  stresses 
iwy  no  longer  be  small  enough  to  bo  ignored.  Referring  to  Fig.  (VII-4), 
the  strains  €  and  e00  are  given  by  the  relations 
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Substituting  Eq.  (VII-7)-(VII-9)  into  Eq.  (VII-46)~(VII-47),  expressions 
for  the  strains  in  terms  of  the  mode  shapes  and  generalized  coordinates 
can  be  obtained.  The  shell  stress  <Jqq  can  then  be  evaluated  from 
Eq.  (VII-43). 
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c)  n  =  0  Partial  Filling  0  <  y  <  1  (See  Figure  VII-5) 


The  major  part  of  the  response  due  to  n  =  0  term3  is  due  to  the  mode 
with  the  lowest  frequency,  j  =  1;  similar  to  the  situation  for  n  /  0. 
Therefore  only  the  one  mode  of  frequency  cr»o,  is  used.  The  radial  displace¬ 
ment  in  this  mode  is 

w  =  qc(t)  (VII-U8) 

while  u(z,  6,  t)  =>  v(z,  9,  t)  =■  0.  The  frequency  and  the  virtual  mass  of 
the  fluid  in  the  partially  filled  tank  are  given  by  Eq.  (V-48),  (V-49), 
and(IV-89),  respectively: 


(VII-52) 


The  generalized  force,  Qo(t),  may  be  written  as  the  sum  of  two  terms: 

Q0(t)  =  Qol  (t)  +  Qo2(t)  *  (VII-53) 

where  QQl(t)  is  the  generalized  force  produced  by  the  radial  blast  pres¬ 
sure  on  the  tank  walls  and  Q^t)  is  the  generalized  force  produced  by 
the  vertical  blast  pressures  acting  on  the  surface  of  the  tank  filling. 

The  force  QQl(t)  is  given  by 

2*  L 

Q0i(t)  =  J  J  He,  Z,  t)  adedz  (VII-54) 

Referring  to  Section  (Vll-a)  and  proceeding  as  in  the  case  for  n  /  ©,  the 

integration  of  Eq.  (VII-51^)  with  respect  to  z  can  be  performed  at  each 

of  the  sixteen  vertical  generators  on  the  shell,  i  =  0,  1,  2,  .  .  .  l6 

as  shown  in  Figure  (VII-l). 

Denoting  the  integral  for  each  i  as  fQi, 

L 

foi  *  /  Pi(z'  t)dz  (VII-55) 

and  letting  P^  be  the  pressure  at  the  various  heights  h  of  the  tank 
(h  =  1,  2,  3, 5,6)  in  the  segment  S  containing  the  point  i,  the  function 
fQi  is  easily  computed  in  terms  of  the  P^.  For  example,  for  the  case  of 
a  9/10  full  tank,  y  =>  0.9/  f  j^ecomes 

foi  "  -3“t  [pu  *  f12  +  pi3  +  •333(P1'»  •  P16>]  (7II'56) 

(6)  For- the  case  of  a  5/10  full  tank,  the  formula  for  f  ^  becomes 

fol  *  -3aS  j^ii  *  pi2  +  P13  +  -333(PiU  -  Pl6)  -  1-333PJ) 
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Equation  (VII-5H)  may  be  written  in  the  form 


Ccl(t)  -  J  V’ 


(VII-57) 


where  f  (0,  t)  is  a  continuous  function  of  the  angle  0  .  The  function 
f  (0,  t)  can  he  expressed  as  a  Fourier  cosine  series 


fo(0>  t}  =  aoo(t)  +  I  aoj(t)cosj0 

j-1 


(VII-58) 


where  the  coefficients  aQ0  and  a  are  determined  so  that  fQ(0,  t)  =  fQl 
when  0  =  0  (i  =  1,  2,  . ...16).  It  will  be  shown  however,  that  only  the 
first  coefficient  aoQ(t)  need  be  computed  in  order  to  determine  Qoi(t). 
Substituting  Eq.  (VII-58)  into  Eq.  (VII-57)  and  integrating,  it  is  seen  that 
only  the  coefficient  aQ0  appears  and  has  to  be  computed.  The  generalized 


force  Qol(l)  ie 


Q  (t)  =  2*a  aQo(t) 


(vii-59) 


where  the  coefficient  aQo  is  obtained  for  n  =  0  from  Eq.  (VII-19). 

The  generalized  force  QQ2(t)  ^ue  to  the  vertical  pressure  acting  on 
the  fluid  surface  at  z  =  7L  is  given  by 
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Qo2(t)  3  ’  £  £  p(r>  9>  t)  ^Q(r)rdrd0 


(VII-60) 


where ^Q(r)  is  the  space  dependent  part  of  ^Q(r,  t),  which  is  the  displace¬ 
ment  of  the  fluid  surface  at  z  *  7L  .  From  Eq.  (lV-91)~(lV-93) : 
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Substituting  Eq.  (VII-6l)  into  Eq.  (VII-60),  QQg(t)  becomes 
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Qo2(t)  =  -  £  £  ?(r>  Q>  ^  IT  +  Y,  PoiJo  rdrd0.(VXI-63) 

w  v,.  u.  1=1 


Referring  again  to  Fig.  (VII-2)  and  (VIX-3),  the  fluid  surface  is 
divided  into  the  subareas  as  in  the  case  n  ^  0.  In  each  subarea, 
the  pressure  is  considered  constant  and  the  integration  of  Eq.  (VII-63) 
can  be  performed  analytically  in  each  subarea.  The  egression  for  QQg(t) 
becomes  therefore 
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The  contribution  to  Qog(t)  from  the  vertical  pressures  acting  on  the 
subarea  S.  in  Eq.  (VII-64)  can  be  vritten  in  the  form 
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where  the  limits  of  the  integration  are  determined  by  the  parti cular  sub- 
area  8,  .  Interchanging  the  integration  with  respect  to  r  and  the 
summation  over  i,  and  defining  the  quantity 
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the  generalized  force  Qq2 


becomes 
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An  analytical  expression  for  IQ^  is  given  in  Table  (VII-l). 


The  generalized  force  QQg(ty  is  given  by 
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which,  upon  substitution  of  Eq.  (VII-67)  becomes 
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As  in  the  case  n  ^  0,  a  sufficient  number  of  fluid  modes  must  be  considered 
commensurate  vith  the  required  accuracy.  For  the  type  of  study  considered 
the  first  three  terms ^of  the  series  of  fluid  modes  suffice  and  the  upper 
limit  in  the  summation  over  i  in  Eq.  (VII-69)  becomes  two. 

The  generalized  force  Qq( t)  can  be  evaluated  from  Eq.  (VII-53)  and 
the  generalized  coordinate  ^(t)  can  be  computed  by  a  numerical  integra¬ 
tion  of  Eq.  (VII-51).  Proceeding  as  in  the  case  where  n  /  0,  and  using 
the  initial  conditions  for  a  Bystem  starting  from  rest,  the  recurrence 
formula  for  qQ(t+k)  becomes 
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X 77  The  constant  displacement  term  is  considered  as  the  lowest  fluid  mode. 
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2«aLm,  1  + 


and  "k"  is  the  interval  of  the  time  steps. 

The  integration  is  started,  using  the  formulas 
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The  generalized  coordinate  ^(t)  gives  the  response  of  the  shell  in 
the  mode  n  =  0.  The  corresponding  hoop  strain,  e00  is  given  by 
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and  the  hoop  stress,  Oqq  becomes 
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a)  n  f  0  Three  Constant  Approximation  -  Partial  Filling 
0  <  7  <  1.  (See  Fig.  VII-5) 


Substituting  the  egression  of  Eq.  (V-3*0~(V-3'>)  f°r  the  mode  shapes 
into  Eq.  (VTI-4)~(VII-6)  for  the  response,  the  displacements  u,  v,  w 
of  a  point  on  the  shell  become : 
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v(z,  0,  t)  -  vo(z)qo(t)  +  Y  cn  f  CO8n0  (VII-77) 

n=l 

Each  mode  is  again  normalized  to  the  total  mass  of  the  empty  tank,  « 
2rtaL  m^  and  the  normalization  coefficients  are  given  by  Eq.  (V-37)* 
The  generalized  coordinates  ^(t)  satisfy  Lagranges  equations, 

Eq.  (VII-10),  vhere  u>n  is  computed  from  the  lowest  root  of  the  frequency 
determinant,  Eq.  (V-33) •  The  components  and  Q^2  of  the  generalized 
force  C^(t)  are  obtained  in  the  same  manner  as  shown  in  Section  (Vll-b). 
For  the  present  case,  the  integral  fni  for  each  point  becomes 
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(VII-78) 


Using  Eq.  (VII -78),  egressions  giving  f^  in  terms  of  the  pressures 
Pih  are  obtained  for  the  cases  of  9/10  full  and  1/2  full  tanks: 

fni  -  V*2  [>a  +  -135Pi2  +  -225P13  +  •095(Pll.  -  Pi6>_  • 

9/10  full  (VII-79) 
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1/2  full  (VII-80) 


Proceeding  exactly  as  in  Section  (Vll-b),  the  generalized  force  component 
Q^(t)  is  obtained  from  Eq.  (VII-24). 

The  generalized  force  q^2(t)  is  given  by  Eq.  (VII-25)  where,  using 
Eq.  (lV-6la)-(lV-62), 
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and  3ni  is  given  by  Eq.  (IV-6lb).  Proceeding  as  in  Section  (VTI-b),  the 
generalized  force  component  ifi  obtained  from  Eq.  (VII-31)*  The 

numerical  integration  formulas,  Eq.  (VII -37) -(VII -39)  can  then  be  used  to 
evaluate  the  generalized  coordinate  q^t) . 

Once  the  coordinates  q^t)  have  been  evaluated,  the  response  of  the 
shell  is  obtained  from  Eq.  (VII-7)-(VII-9) •  Although  the  three  constant 
approximation  may  be  used  to  obtain  an  estimate  of  the  response  of  the  shell, 
it  does  not  give  accurate  results  for  the  strains  and  stresses  in  the  shell. 
It  must  be  emphasized  that  when  shell  strains  and  stresses  are  required. 


the  five  constant  approximation  must  be  used. 
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e)  n  -■£  0  Empty  Shell  -  Five  Constant  and  Three  Constant  Approximations 


In  Section  (VH-b)  and  (Vll-d),  methods  for  the  determination  of  the 
transient  response  of  partially  filled  fuel  tanks  to  blast  loadings  are 
developed.  The  methods  presented  in  these  sections  are  general  and  can 
be  applied  to  the  case  of  empty  tanks,  provided  that  the  correct  normalized 
modes  and  frequencies  of  empty  tanks  j  Sec.  (I-a)  and  (i-e)  J  are  used  in 
the  analysis.  Moreover,  the  generalized  force  component,  Qfa2(t);  must  be 
set  equal  to  zero  for  the  case  of  an  empty  tank.  With  these  exceptions, 
the  procedures  of  the  previous  sections  can  be  used  in  determining  the 
transient  response  of  empty  cylindrical  fuel  tanks  to  blasts. 

Specific  formulas  in  terms  of  the  wall  pressures  cannot  be  given 
at  this  time,  since  the  Armour  report  (See  footnote  1,  Section  VIl) 
contains  only  pressure  information  for  9/10  full  and  1/2  full  tanks. 
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f)  n  =  0,  Empty  Shell. 

The  response  of  empty  shells  for  the  special  case  n  =  0  haB  been 

discussed  in  Section  (i-e)  and  (Vll-a).  To  obtain  the  contribution  of 

the  mode  n  =■  0,  the  empty  shell  is  considered  to  act  as  a  series  of 
(7) 

separate  rings'  }  each  ring  is  subjected  to  the  radial  pressure 
acting  over  a  circumferential  strip  of  shell  equal  in  depth  to  the  depth 
of  the  vertical  subarea  "ih"  (_See  Fig.  (VII-l)J  .  The  response  of 
each  ring, 


woh  ■ 


(VII-82) 


satisfies  Lagrange s  equations 


%h  +  ^oh  %h 


M 


(VII-83) 


where  <0^,  the  frequency  of  each  ring,  is  given  by 
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and  Mq  7  the  generalized  mass,  is  given  by 
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The  generalized  force  on  each  ring  is  denoted  by  Q^t)  and  can  be 

obtained  in  a  manner  similar  to  that  of  Section  (VII-c).  Letting  P^(t) 

be  the  radial  blast  pressure  at  the  various  points  i  on  the  particular 

ring  defined  by  the  index  h,  the  generalized  force  becomes 

2n 
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(VII-86) 


TT5  Each  ring  is  identified  by  an  index  "h",  corresponding  to  a  particular 
subarea  "ih",  (jSee  Fig. (VII-l)J 


where 


8 

ph(S,  t)  -  aQo(t)  +  £  aoJ(t)cosj0'  (VII-87) 

j=l 

The  coefficients  &0y  J  =  0-8,  are  determined  so  that  t)  =  P^^l) 

when  fl  s  fl  (i  e  l,  2,.  .  .  .  16).  Substituting  Eq.  (VII-87)  into 
Eq.  (VII-86)  and  integrating,  it  is  seen  that  only  the  coefficient 
aoo(t)  appears  and  has  to  be  cornput  <1.  The  generalized  force  (^(t) 
is 

Qtoh(t)  =  2«a  aoQ(t)  (VII-88) 

where  the  coefficient  a  (t)  is  given  by 

oo 

15 

aco(t)  ”  35  X  V^S5  (VII-89) 

i=0 

Proceeding  as  in  Section  (VII-c),  Eq.  (VII-83)  can  be  integrated 
using  Eq.  (VII-70)  and  (VII-72)  where 

a  (t) 

F  (t)  =  -°£ -  .  (vii-90) 

i 

The  response  w^t)  is  assumed  to  act  at  the  center  of  each  ring  and  the 
corresponding  hoop  strain  and  hoop  stress  for  each  ring  can  be  evaluated 
using  Eq.  (VII-73)“(VII-7^)  respectively. 
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Appendix  to  Section  V.tl  -  Determination  of  the  Fourier  Series  Coefficients  a 


The  continuous  function  f^{0,  t)  given  by  Eq.  (VII-18), 

8 

fn(0,  t)  -  anQ(t)  +  an<j(t)  cosjS 


(A-l) 


is  required  for  the  determination  of  the  generalized  force  component  Q  ..(t) 

The  coefficients  a  .  are  chosen  so  that  f  (0,  t)  =  f,  when  0  «*  0. 

nj  n  ni  l 

(i  =  l-l6).  In  this  section,  the  expressions  for  the  coefficients  a  in 
terms  of  the  quantities  f  will  bo  developed. 

Consider  the  Fourier  Series 


f  (0)  =  a  +  )  a  .  cosj0 
n'  '  no  nj 


(A-2) 


where  the  coefficients  a  are  given  by  the  expressions 
2n  ^ 


J  fn(e)  -osje  de 

o _ 

r2rt  2 

/  cos  J0d0 


(A-3) 


Let  the  value  of  the  function  fn(0)  he  known  at  N  equally  separated 

points  i,  1  =  l-l6,  where  each  point  i  is  located  by  the  angle  0 (i  ■=>  l, 

2,  3.  .  .N).  The  coefficients  a^  of  the  series,  Eq.  (A-2),  are  to  be 

evaluated  so  that  f  (0)  =  f  (0  )  =  f  .  at  each  point  i.  Calling  the 

n  n  1  ru 


interval  of  separation 


(A-4) 


the  coefficient  of  Eq.  (A-3)  can  be  written  as  a  ratio  of  the  summation 

,  N 

over  i .  t— > 


Setting  9i  =  — ! -  ,  1  -  0,  1,  2, 


Becomes:  N-l 


V"  *  /2nix  ,  2ni 

l  f„i(— )co8j 


i=0 


nj  N-l 


I 

i=0 


2  2rti 
C°B  J  ~ 


.N-l,  for  convenience,  Eq.  (A-5) 


(A-6) 


Using  the  summation  formula 
N-l 


Z2  ,  N  sina(2N-l)  +  since, 
cos  ai  =  p  +  - 1 - 1 - 

1-0  ^oina 


and  setting  a  =  j,  the  following  results  are  obtained: 


(A-7) 


N-l 


4/^4/N  V  2  2n  . .  N 
o  f  0,  J  +  t,  :  )  cos  —  Ji  =  - 

d  l_j  ii  2 

i=0 


(A-8) 


J  =  0 


J  = 


N-l 

I 

i=0 

N-l 

I 

i=0 


cos  oi  =  N 


cos  >ri  =  N 


(A-9) 


(A-10) 


Substituting  Eq.  (A-8)-(A-10)  into  Eq.  (A-6),  the  coefficients  a^  become 


nj 


2  V  ^  /2«1n  i  2nl 

-  5  l  fnl(— >  COOJ 


N 


i=0 

N-l 


no 


N  ni^N  ' 


(A-ll) 


(A- 12) 


i=0 

N-l 


l  V  f  /2nix 
an,N  N  (_>  ni'  N  ' 
2  i=0 


COBjti 


(A-13) 


For  the  particular  case  where  !?  -  lo,  Eqs.  (VII-19)-(VII-2l)  are 
obtained.  It  should  be  noted  that  when  a  harmonic  analysis  passing  a 
Fourier  Series  through  N  points  is  made,  the  series  terminates  after 
the  term  (|  +  l). 
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n  =  6 


a  4r 

r  J6<4->  dr 
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ni 


144a2- a2  r2  a  r 

+  (-s—Si-)  ^(4-) 


tt  .T 

rJ6(4-)ar  . 
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ni 
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2^6a  +  anir 

a  .  r  a 
ni 


00 
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(*)  For  n  =  odd  integers,  the  integral  is  obtained  as  an  infinite  series 


in  ^2n+l^ani  r/a)  a  num^er  'terms  commensurate  with  the  required 
accuracy  must  be  taken.  Moreover,  the  integral  between  limits  r^  and  r2  is 


computed  from: 


«  2  r  2  r‘l 

/  f(r)  dr  a  /  f(r)dr  -  /  f(r)dr. 
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FIG.M-2  VERTICAL  PRESSURE  ON 
FLUID  SURFACE  ARRANGEMENT 
AND  NOMENCLATURE 


FIG.M-3  TYPICAL  SECTOR 

FOR  VERTICAL  PRESSURES 
IN  FLUID  SURFACE 


FIG.SE -4  NOMENCLATURE 
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VIII  CONCLUSION  (Comment  on  the  Purposes  and  Results  of  the  Mode  Analysts) 

The  logical  first  step  in  the  analytical  study  of  a  new  problem  concern¬ 
ing  the  ability  of  a  structure  to  withstand  static  or  dynamic  loads,  is 
the  formulation  of  a  theory  giving  the  elastic  (usually  small  deflection) 
response  of  the  structure.  Such  a  theory  will-in  general-not  be  valid  up 
to  the  stage  where  the  structure  collapses  or  loses  its  usefulness  otherwise] 
instead,  the  purpose  of  this  simplest  theory  is  to  find  out  which  of  the 
many  possible  refinements,  large  deflections,  plastic  behavior,  etc.,  ought 
to  be  included  in  an  analysis  to  make  it  valid  up  to  capacity  loads. 

An  obvious  exaaple  for  the  situation  outlined  is  the  capacity  of  beams 
to  carry  static  loads i  the  conventional  (elastic)  theory  ceases  to  be 
valid  when  the  yield  point  is  reached  somewhere  in  the  cross  section.  To 
obtain  a  theory  valid  up  to  capacity  loads  plastic  effects  must  be  included. 
It  is,  however,  not  always  a  plastic  theory  which  is  required,  as  may  be 
seen  from  the  example  of  a  slender  column.  In  this  well-known  case,  elastic 
second  order  effects  must  be  included  to  be  able  to  describe  the  buckling 
behavior. 

Similarly  in  the  present  study  of  the  ability  of  fuel  tanks  to  resist 
air  blast,  the  mode  analysis  cannot  be  expected  to  be  valid  up  to  the  end 
of  the  usefulness  of  the  structure,  because  the  a-priori  unknown  lethal 
damage,  possibly  tearing  of  the  plates,  buckling,  sloshing  of  fluid,  nec¬ 
essarily  will  be  preceded  by  a  phase  requiring  one  or  more  refinements 
beyond  an  elastic  small  deflection  theory.  The  mode  analysis  presented 
can  however  be  used— and  has  been  formulated  with  considerable  effort 
expressly  for  the  purpose  of  determining  which  physical  effects  will 
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occur  requiring  a  new  theory.  The  sole  objective  of  the  mode  analysis  is 

therefore  to  find  the  end  of  its  own  usefulness,  occur^.ng  at  a  certain 

time  after  a  blast  wave  hits  the  tank,  and  disclosing  which  additional 

physical  effects  must  be  Included  in  a  subsequent  stage. 

The  stresses  in  a  typical  100  ft.  diameter  tank  (steel,  9/10  full) 

were  determined  by  means  of  this  mode  analysis,  and  are  presented  as 

sample  computation.^  In  addition,  less  detailed  determinations  of 

stresses  were  made  on  full,  half  full  and  empty  full-scale  protected  and 

unprotected  tanks,  and  on  scaled  model  tanks  ranging  from  l6"-72"  in 

diameter.  These  computations  indicate  in  all  cases  that  the  mode  analy- 

(2) 

sis  ceases  to  be  valid  after  a  quite  short  time, '  the  reason  being  that 

the  tank  wall  on  the  side  facing  the  wave  lifts  from  the  support. 

The  mode  analysis  assumes  that  the  shell  is  cantilevered  from  the 

base  and  that  the  shell  is  prevented  from  vertical  (and  horizontal) 

motion,  an  assumption  which  is  initially  correct.  However,  the  support 

of  the  shell  iB  such  that  large  downward  reactions  can  be  carried,  while 

uplift  of  the  shell  is  prevented  only  by  the  dead  weight  of  the  shell 

(and  possibly  by  nominal  bolts  of  small  total  strength).  E.g.,  in  a 

steel  shell  of  100  ft.  diameter  and  40  ft.  height,  the  vertical  compressive 

stress  is  only  340  lbs. /in  ,  such  that  any  vertical  tensile  stress 
2 

o  >  340  lbs. /in  ,  induced  by  the  blast  will  cause  uplift.  At  load  in- 
tensities  of  interest,  vertical  tensile  stresses  of  such  small  magnitude 
develop  very  quickly  and  the  validity  of  the  movie  analysis  ends  due  to 

(1)  See  Part  (III)  of  the  Final  Report 

(2)  Short  in  comparison  to  the  decay  constant  of  the  shockwave. 


uplift  for  a  100  ft.  diameter  tank  after  only  a  few  millicccondrMlong  before 
plasticity  or  large  deflections  effects  can  develop. 

'Die  theoretical  conclusions  that  uplift  will  develop  has  since  been 
verified  by  teats  on  mercury  filled  16  inch  diameter  models  in  the  shock 
tube,  and  by  teats  on  wuter  filled,  full,  half -full  and  empty  models  of 
10  ft.  diameter.  The  testa  have  demonstrated  (See  Fig.  VIII-l)  clearly, 
that  large  uplift,  exceeding  25  per  cent  of  the  height  of  the  model  does 
occur  on  models  prior  to  lethal  damage.  This  fact  leads  to  the  conjecture 
that  appreciable  uplift  may  also  occur  on  full  scale  tanks. 

From  the  instant  of  time  when  the  mode  analysis  becomes  invalid,  a 
new  type  analysis  is  therefore  required]  this  analysis  is  complicated 
by  the  fact  that  a  part  of  the  shell }  (Fig.  VIII-2);  over  some  unknown  and 
possibly  varying  angle  2a  will  not  be  vertically  supported,  while  the 
remainder  of  the  shell,  wherever  the  reactions  are  downward  is  still 
supported  by  the  base.  Analytical  treatments  for  this  situation  are 
developed  in  Part  II  of  this  Report.  These  treatments  are  also  made  in 
stages,  because  the  first  theory  allowing  for  uplift  only  requires  addi¬ 
tional  refinements  in  order  to  apply  up  to  lethal  domage. 

In  view  of  this  fact  that  tentative  computations  and  model  tests 
indicate  that  the  mode  analysis  always  ceases  to  be  valid  long  before 
lethal  damage  develops,  it  appears  that  the  mode  analysis  need  not  be  re¬ 
peated  at  all  in  a  study  of  effects  of  various  blast  waves  on  tanks  of 
(3) 

different  sizes,  ,  but  that  such  a  study  requires  only  the  application 
of  the  methods  developed  in  port  II.  It  is  noted  that  the  work  on  the 


(3)  Unless  some  unconventional  type  tank  is  studied. 


mode  analysis  recorded  in  Port  I  not  only  was  required  to  recognize  that 
tanks  will  lift  prior  to  failure,  but  that  a  large  part  of  the  derivations 
are  used  again  in  the  theories  developed  in  Part  II. 


The  problem  of  the  static  buckling  of  an  empty  cylindrical  fuel  tank 

subjected  to  a  unform  radial  pressure  producing  a  uniform  compressive 

hoop  stress  and  strain  €s  considered  in  this  Appendix.  ^  The  fuel 

tank  is  of  height  L  and  radius  a,  and  is  free  at  the  top  and  simply 

supported  at  the  base.  The  uniform  buckling  pressure,  p^,  is  obtained 

by  use  of  the  Rayleigh-Ritz  method  which  leads  to  a  direct  solution  of 

ths  extremum  problem  arising  from  the  theorem  of  Stationary  Potential 
(o\ 

Energy.  '  This  theorem  is  given  by  the  following  relation: 

II  ■  ?  +  W  »  stationary  •  (B-l) 

4 

where  the  stationary  value  of  the  total  potential  energy  U  must  be  a 
minimum  for  stable  equilibrium.  The  quantity  V  represents  the  strain 
energy  of  the  structure  and  W  is  the  potential  energy  of  the  external 
loads. 

The  general  expression  for  the  strain  energy  V  in  terms  of  the  shell 

4' 

displacements,  u(z,  ©),  v(z,  9)  and  v(z,  9),  is  given  by  Eq.  (I-10)-(I-12) 
of  the  main  paper  which  are  repeated  here  for  convenience: 

V  -  \  +  V2  (B-2) 

(1)  This  problem  was  considered  because  Of  the  possibility  that  tanks 
might  buckle  due  to  the  hoop  stress  produced  by  the  blast  loads.  It  was 
Ofound  subsequently,  that  in  the  case  of  interest,  uplift  of  the  tanks 
occurs  prior  to  any  possibility  of  buckling  so  that  the  analysis  reported 
here  was  not  utilized. 

(2)  "Buckling  Strength  of  Metal  Structures’’  by  Friedrich  Bleich,  Edited 
by  Hans  n.  Bleich,  McGraw-Hill  Book  Company,  1952,  New  York,  Pg  TO  ff. 
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where  the  first  term  to  llie  membrane  attain  energy : 

2n  L 


Vl" 
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2(1*^) 


h 
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2  2 
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o  o 


u,;  +  (ve  -  mT  +  2u\)u z(v0  ~  w) 

+  (^r)(up  +  avz)2J  dzd<? 


(B-3) 


and  the  second  term  Vg  represents  the  strain  energy  of  bending^ and  coupling 
terms  between  the  membrane  and  bending  strains: 
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.  w  n,  +  2a  vw  (vfln  +  v„) 
z  z0'  zzv  00 

dzd0  .  (B-4) 


Tli-  expression  for  W,  the  potential  energy  of  the  external  uniform 

(3) 


pressure  pn,  is  given  by 

il 


L  2n 


2  n 
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/  |y  -  4  -  4  -  2awuz  _ 

Pna  r  1 

as  az  -  -g-  J  vu0 
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(B-5) 


Using  the  Rayleigh-Ritz  method,  an  appropriate  set,  of  coordinate 
functions  with  arbitrary  coefficients,  which  satisfy  the  boundary  condi¬ 
tions  on  the  deflections  u,  v,  and  w  of  the  shell  is  assumed-  Substitut¬ 
ing  these  expressions  into  Eq.  (B-l) -(B-5),  the  total  potential  energy  U 
is  obtained  as  a  function  of  the  arbitrary  coefficients.  For  static 
equilibrium,  the  condition  that  U  be  a  minimum  leads  to  on  ordinary  maximum- 
minimum  problem  on  the  arbitrary  constants  and  hence  to  a  set  of  homogeneous 


(3)  "Eynamic  Response  of  Cylindrical  Tanks"  by  F,  L«  DiMaggio,  Armed 
Forces  Special  Weapons  Project,  Contract  DA~29”044~XZ~54,  AFSWP  No,  1075, 
May  1956,  Pg.  21,  Eq.  (6). 


linear  algebraic  equations  on  these  constants.  For  non  zero  solutions, 
the  determinant  of  the  set  of  homogeneous  equations  must  be  set  equal 
to  zero,  thus  leading  to  an  equation  on  the  critical  pressure  p 
corresponding  to  a  buckling  mode  defined  by  the  number  "n"  of 
circumferential  waves  in  the  shell  displacements.  Buckling  determinants 
will  be  developed  for  both  the  three  constant  and  the  five  constant 
approximations  considered  in  Sections  (i)-(VIIl). 
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(u)  Three  Constant  Approximation 


Let  the  di splacemento  of  the  ehell  (Fig.  I-l)  be  given  by  the 
coordinate  functions 


u(z, 

e)  =  u  - 

'  a 

cos  n0 

(B-6) 

v(z. 

0)  =  V  - 
a 

sin  n0 

(B-7) 

w(z, 

0)  =  W  - 
'  a 

cos  n0 

(b-8) 

where  U,  V  and  W  are  undetermined  constant  coefficients.  Substituting 
into  Eq.  (B-2)-(l3-5)>  the  total  strain  energy  V  is  given  by 
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and  the  potential  energy  of  the  external  pressure  W  becomes: 


lillllj2  w2  -  U2  -  |WU  +  njvuj 


(B-10) 


The  condition  of  minimum  total  potential  energy 


(B-ll) 


where  successively  takes  the  values  U,  V,  W  leads  to  the  following  set 
of  three  homogeneous  linear  equations  in  the  unknowns  U,  V  and  W: 
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where 


Eh 


(B-15) 


v.  “ 

and  k  =  - _  .  Non  zero  solutions  of  Eq.  (B-12)-(B-l4)  require  that  the 

12a 

determinant  of  the  system  vanishes,  thus  leading  to  the  determenantal  buckl¬ 
ing  equation  shown  on  Page  165 . 

The  lowest  root  M  _  of  Eq.  (B-l6)  defines  the  critical  buckling  pres- 
ild 

sure  p  „  in  the  mode  denoted  by  n.  The  buckling  stress  a  „  can  then 
UB  Hd 

be  computed  from  the  relation 


nB 


a 

PnB  h 


(B-17) 


It  should  be  noted  that  for  very  thin  shells  in  which  ~  is  very  small, 

cl 

the  constant  k  can  be  set  equal  to  zero  in  Eq.  (B-l6)  as  in  the  vibration 
problem  of  Section  (l). 


Figure  (l3-l)  shows  a  curve  of  the  buckling  load  number  plotted 
against  the  mode  number  n  for  an  unprotected  steel  tank  (no  concrete 
shielding).  For  the  lower  modes  l  <  n  <  7>  the  shell  acts  as  a  membrane 
while  for  the  higher  mode 3  n  >  18,  the  bending  strains  predominate  and 
control  the  critical  buckling  loads. 

As  an  example,  the  following  table  shows  the  static  hoop  stresses 
for  n  =  2,  6,  16  at  which  a  steel  tank  of  100  ft.  diameter,  1*0  ft.  height 
and  1/2  inch  wall  thickness  would  buckle.  The  smallest  hoop  stress  is 
required  for  n  =  16. 

n  2  6  l6 

o^psi)  210,000  17,81*0  720 

Figure  (B-2)  shows  the  similar  curve  for  a  protected  steel  tank 
(18  inch  concrete  protection).  For  protected  steel  tanks,  the  buckling 
pressures  png  are  of  course  very  much  higher;  in  addition  the  minimum 
occurs  for  a  much  smellier  value  of  n. 


(b)  Five  Constant  Approximation. 


Let  the  displacements  of  the  shell  be  given  by  the  more  general 


coordinate  functions 


u(z,0) 


cos  nd 


(B-18) 


v(z,  0 ) 


V  -  sin  n0 
a 


(B-19) 


w(z,  0) 


1  +  «  <t  -  K> 


cos  nd 


(B-20) 


vhere  U,  V,  W,  X  and  Y  are  undetermined  constant  coefficients. 

Proceeding  exactly  as  in  Part  (a)  of  this  Appendix,  the  determenantal 
buckling  equation  shown  on  Page  168  is  obtained. 

The  lowest  root  MQg  of  Eq.  (B-21)  defines  the  critical  buckling 
pressure  p^  in  the  mode  n. 

In  sample  computations,  it  was  found  that  the  buckling  loads  are  well 
approximated  by  the  three- constant  approximations  and  that  it  is  unnecessary 
to  use  the  fifth  order  determinant,  Eq.  (B-21).  This  can  be  seen  from 
Fig.  (B-2). 
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(c)  dynamic  Duckling 


The  term  "dynamic  buckling"  is  only  another  way  of  saying  that 
the  structure  will  show  a  transient  response.  This  transient  response 
can  be  found  again  from  the  appropriate  equations  of  motion  using  the 
modes  of  free  vibration  as  generalized  coordinates.  In  the  presence  of 
a  hoop  stress  a  which  is  a  substantial  percentage  of  the  buckling  stress 
anB*  °r  excee<i8  °rlg^  Eq.  (V.'.'I-10)  the  equation  of  motion  for  the 


generalized  coordinate  becomes 


•  •  2  /  ,  Ox  Sn 

%  +  “n  (  1  '  JT  >  %  =  =- 
nB 


(B-22) 


If  o/o  is  small,  such  that  it  can  be  neglected,  the  original  equation 
'  nb 

of  motion,  Eq.  (VII-10)  is  again  obtained.  It  is  noted  that  Eq„  (B-22) 
is  valid  for  both  positive  or  negative  values  of  1  - —  . 
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